GLOBAL SOLUTIONS FOR THE GRAVITY WATER WAVES EQUATION IN 

DIMENSION 3 



P. GERMAIN, N. MASMOUDI AND J. SHATAH 



Abstract. We show existence of global solutions for the gravity water waves equation in dimension 
3, in the case of small data. The proof combines energy estimates, which yield control of related 
norms, with dispersive estimates, which give decay in L°° . To obtain these dispersive estimates, we 
use an analysis in Fourier space; the study of space and time resonances is then the crucial point. 



1. Introduction 

We consider the three-dimensional irrotational water wave problem in the presence of the grav- 
ity g. The velocity of the fluid is denoted by v. We assume that the domain of the fluid is 
given by = {{x^z) = {xi,X2,z) E ]R.^,z < h{x,t)}, where the graph of the function h, i.e. 
S = {{x,h{x,t)),x £ M^} represents the free boundary of the fluid which moves by the normal 
velocity of the fluid. In this setting the Euler equation and the boundary conditions are given by 



(E) 
(BC) 



def 

Dtv = dtv + VyV = -Vp - ges {x, z) G O, 
V • ^; = {x,z) ^VL, 



\dth + v ■V(^^^^){h- z) = xGM^ 
\p\s = 0. 

Since the flow is assumed to be irrotational, the Euler equation can be reduced to an equation 
on the boundary and thus the system of equations (E-BC) reduces to a system defined on S. This 
is achieved by introducing the potential ipj-i where v = VV'w- Denoting the trace of the potential 
on the free boundary by '4'{x^t) = ip-j-i^x , h{x , t) , t) , the system of equations for ip and h are [SS99j 

r dth = G{h)i; 

(WW) <^ dti^ = -gh - ^ |V^p + 2(i+ivfei^) iG{h)i^ + Vh • Vi^y 

( (/i,V)(t = 2) = (/io,Vo). 
where G{h) = y^l + |V/ipAA, J\f being the Dirichlet-Neumann operator associated with Q. 

def , 

Before stating our result we introduce the Calderon operator A = \D\, the complex function 



\2 



U 



h + iA^I'^il), its initial data uq ho + iA^^'^ipQ, and its profile / e^^^^^^u. Then 



Theorem. Let 5 denote a small constant and N a big enough integer. Define the norm 

\\u\\x sup t \\u\\y^^4,oo + t^^ \\u\\hn + \\xf\\]^2 + \\u\\l2 . 
t>2 

There exists an e > such that if the data satisfies ||e**^ ^ uq\\x < s, then there exists a unique 
global solution u of (WW) such that \\u\\x ^ £• 

The estimates leading to the above theorem give easily that the solutions scatter. 

Corollary 1.1. There exists a constant Cq such that: under the conditions of the previous theorem, 
there exists foo G H^-^CoS p ^2(^2-005^^-, ^^^^^ ^^^f j^^^ _^ ^N-NCoS p L^^x^-CoS^^^ 

as t ^ oo. 
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Some comments on the stated results. Our approach in writing this article has not been to 
strive for a minimal A^; rather, we chose to give a proof as concise as possible, by allowing at some 
places losses of derivatives. The fact that the data for (WW) are given at i = 2 does not have a 
deep significance, it is simply a matter of convenience to avoid writing 1/(2 + t) when performing 
estimates, since the decay of 1/i given by the linear part of the equation is not integrable at 0. 
The precise conditions on the constants 6, N, e under which the theorem holds are, as an inspection 
of the proof reveals: 6 less than a universal constant 5o, ^ <5 and e 5. Finally, we notice that 
the condition on uq in the theorem is satisfied if, for instance, uq G W^'^ D n L?'{x'^dx). 

There is an extensive body of literature on local well posedness and energy estimates for (WW) 
in Sobolev spaces starting with the work of Nalimov |NA74j for small data in 2 dimensions (see 
also H. Yoshihara [YQ82j ). The first breakthrough in solving the well posedness for general data 
came in the work of S. J. Wu |WU97t IWU99j who solved the problem in 2 and 3 dimensions. 
There are many other works on local well posedness: we mention the work of W. Craig |CR85j . T. 
J. Beal, T. Hou, and J. Lowengrub [BHL93], M. Ogawa and A. Tani [OT02|. G. Schneider and E. 
C. Wayne |SW02j . D. Lannes [LA05j . D. Ambrose and N. Mamoudi [AM06] [AM07j . and P. Zhang 
and Z. Zhang |ZZ06| . We also mention that for the full system (E) there are well posedness results 
given by D. Christodoulou and H. Lindblad [CLOOj . D. Coutand and S. Shkoller |CS05] . J. Shatah 
and C. Zeng [SZ06j. Recently S. J. Wu proved almost global existence for small data for (WW) in 
two space dimensions |WU09j . 

Our proof of global existence is based on the method of space-time resonances introduced in 
[GMS08) . It is the key notion to understand the long time behavior of solutions to the water waves 
problem. The space time resonance method brings together the normal forms method [SH85] and 
the vector fields method |KL85j . and is illustrated below. 

Space-time resonances. If one considers a nonlinear dispersive equation on M" 

dtu = iL(^-\7^u + au^ + . . . , 
solutions can be represented in Fourier space in terms of Duhamel's formula 

(1.1) fit, = MO + fj e*^^«'^)a/(5, 7?)/(s, i-r,)+ ... d^ds, 

where (p = — i(C) + L{r]) + — r]), /(^) = e~^^^^^^u, and where u denotes the Fourier transform 
of u. The reason for writing the equation in terms of / rather then u is that / is expected to be 
non oscillatory and thus for the quadratic terms the oscillatory behavior is restricted to e^'^'^^^'^h 

• Time resonances correspond to the classical notion of resonances, which is well known from ODE 
theory. Thus quadratic time resonant frequencies are defined by 

=^ = {(e, r/); vie, V)=0} = {(e, r/); L(r/) + L(e - v) = HO}. 

In other words, time resonances correspond to stationary phase in s in (jl.ip . In the absence of time 
resonant frequencies, one can integrate by parts in time to eliminate the quadratic nonlinearity in 
favor of a cubic nonlinearity. The absence of time resonances for sufficiently high degree usually 
indicates that the asymptotic behavior of solutions to the nonlinear equation is similar to the 
asymptotic behavior of linear solutions. 

However, for dispersive equations time resonance can only tell part of the story. The reason is 
that time resonances are based upon plane wave solutions to the linear equation e*(^(^)*~*''>'^\ while 
we are usually interested in spatially localized solutions. For this reason we introduced the notion 
of space resonances. 

• Space resonances only occur in a dispersive PDF setting. The physical phenomenon underlying 
this notion is the following: wave packets corresponding to different frequencies may or not have the 
same group velocity; if they do, these wave packets are called space resonant and they might interact 
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since they are locahzed in the same region of space. If they do not, they are not space resonant and 
their interaction will be very limited since their space-time supports are (asymptotically) disjoint. 
If one considers two linear solutions ui and U2 whose data are wave packets localized in space 
around the origin and in frequency around rj and — rj, respectively, then the solutions ui and U2 
at large time t will be spatially localized around {—VL{r])t) and (— V-L(^ — T])t). Thus quadratic 
spatially resonant frequencies are defined as 

y = {(e,r/); VL(r?) = VL(^ - r?)} = r/); V^(^(e, r?) = 0}; 

in other words, space resonances correspond to stationary phase in tj in (jl.ip . Spatial localization 

can be measured by computing weighted norms, namely ||x^u(0, •) 

or = •) for the solution. 

This leads to the observation that if there are no space resonant frequencies then one can capture 
the spatial localization of the solution by integrating by parts in rj in equation (jl.ip . 
• Space-time resonances correspond to space and time resonances occurring at the same point in 
Fourier space. The space-time resonant set is given by 

It is clear from the above description that the asymptotic behavior of nonlinear dispersive equa- 
tions, of the type considered here, is governed to a large part by the space-time resonant set 
Thus if the set M ^ ^ and the dispersive time decay is weak, then one faces serious difficulties in 
proving that small solutions exist globally and behave asymptotically like linear solutions. How- 
ever if the quadratic nonlinearity vanishes on M then we say that the quadratic term is null and 
try to eliminate it by integration by parts in time or frequency. After dealing with the quadratic 
nonlinearity we move on to study the cubic interactions where we have similar definitions of 
and ^. 

When applying this method to the problem on hand, (WW), one can easily compute that the 
quadratic terms are null; as for the cubic terms, their resonance structure is more subtle, but they 
can be estimated by the same method. However the integration by parts procedure, in addition 
to introducing multipliers with non smooth symbols, introduces new singularities due to the non 
smoothness of the dispersion relation which is given by |^|2. After dealing with these difficulties, 
the problem is reduced to studying the asymptotic behavior of quartic and higher order terms. But 
for quartic terms the dispersion is strong enough to overcome any effect of resonance, thus we are 
able to prove global existence. 

Plan of the article. In order to prove the existence result, we will prove the a priori estimate 
||u||x < 00. The organization of the article is as follows: 

Section\^ contains energy estimates (control of the part of the X norm), using the formalism 
developed in |SZ06j . Most of the rest of the paper is dedicated to controlling the L?'{x'^dx) and 
L°° parts of the X norm. Section is the first step in this direction: the space-time resonant 
structure of the quadratic and cubic terms in the nonlinearity of (VFI^) is studied. In Section^ a 
normal form transform is performed. This leads to distinguishing four kinds of terms: quadratic, 
cubic weakly resonant, cubic strongly resonant, and quartic and higher-order, or remainder terms. 
Quadratic terms are treated in Section [5l weakly resonant cubic terms are treated in Section O 
strongly resonant cubic terms are treated in Section and quartic and higher-order terms are 
treated in Section\^ Finally, scattering (Corollarv ll.ip is proved in Section\^ 

The appendices contain more technical developments needed in the main body. Some classical 
harmonic analysis results are recalled in appendix H] and results on standard pseudo-product oper- 
ators are given in appendixWl Particular classes of bilinear and trilinear pseudo-product operators 
are defined and studied in appendices [0 and O respectively. Estimates on the Dirichlet-Neumann 
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operator are proved in appendix and finally appendix [71 gives bounds on the quartic and higher 
order terms in the nonlinearity of (WW). 

Notation. Throughout this manuscript we will use the following notations. 

A < B oi B > A means that that A < CqB for a constant Cq. 
A-^ B means that A< B and B < A. 

A « B or B >> A means that A < eqB for a small enough constant eq. 

We denote by c (respectively C) small enough (respectively big enough) constants whose values 
may change from one line to the other. We use standard function spaces; their definitions are (see 

appendix lAl for the definition of the Littlewood-Paley operators Pj) 

1/9 



def 



Homogeneous Besov spaces Bp ^: 

Homogeneous Holder spaces C°', with < a < 1: equivalent to B, 

def 



a 

oo,oo ■ 



l^<0/||, 



Inhomogeneous Besov spaces Bp ^: ||/||_b| 

Inhomogeneous Sobolev spaces PV^^'P: H/lliy^.p = II/II-Lp 
Inhomogeneous Sobolev spaces H^: equivalent to W^'"^ or -62^2- 



1/g 



We also make use of the standard embedding relation between these spaces, and first of all of the 
Sobolev embedding theorem: \il <p < q < 00 and /c > | — |, then < H/Hi^/fc,?. 

The energy estimate will be reproduced along the lines of |SZ06] . thus in deriving the energy 
estimates we will adhere to the same notation as in [SZ06J. 

A*: the adjoint operator of an operator; Ai ■ A2 = trace(Ai(^2)*)j for two operators. 

D and d: differentiation with respect to spatial variables. 

Vx- the directional directive in the direction X. 

Dt = dt + v^d^i the material derivative along the particle path. 

S = {{x, h{x)) G M^; /i : — > M} a surface in given as the graph of h. 

= {(x, z) G M^; z < h{x)} the region in bounded by S. 
N: the outward unit normal vector to = dil.. 

_L and T: the normal and the tangential components to S of the relevant quantities. 

T^w = V^, for any w tangent S; T> the covariant differentiation on S* C M^. 

H: the second fundamental form of 5, Ii{w) = V^N; U{X, Y) = U{X) ■ Y. 

k: the mean curvature of S, i.e. k = trace H. 

f-H- the harmonic extension of / defined on S into 0. 

■^if) = ^ nIh : 5* ^ M: the Dirichlet-Neumann operator. 

A5 = traceP^: the Beltrami-Lapalace operator on S. 

(-A)-^ is the inverse Laplacian on Q, with zero Dirichlet data. 

The decay and weighted estimates will be carried out using harmonic analysis techniques, where 
we employ standard notations and denote by: 

/ or T f: the Fourier transform of a function /. 

m{D)f = J^~^in(^){J^f ){£,)'■ the Fourier multiplier with symbol m. In particular A = \D\. 

, Xfi)- a smooth scalar or vector-valued function in all its arguments (xi, . . . , x„). 
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A{xi, . . . , Xn)[y]- a function which is hnear in y and smooth in (xi, . . . , 

A{xi, . . . Xn)[y, z\. a function which is bihnear in (y, z) and smooth in (xi, . . . , x^)- 

2. Local existence and Energy Estimates 

There are a variety of methods that have been developed to prove local well posedness and to 
obtain energy estimates for (E-BC) in Sobolev spaces. For example in [SZ09j local well posedness 
was established in K'^ for m > 5/2. The success of all of these methods hinges upon finding the 
appropriate cancellations needed to isolate the highest order derivative term with the appropriate 
Rayleigh- Taylor sign condition. Here we follow the method developed in [SZ06j to obtain higher 
energy estimates. The idea is the following: since the system is reduced to the boundary we apply 
the surface Laplacian to (BC) instead of partial derivatives. This leads to an evolution equation 
for the mean curvature k. 

Conserved energy. We start by deriving the conserved energy of the (E-BC) system. Observe 
that the Euler system in the presence of gravity can be written as 

dtv + V ■ Vv = —Vp — ges —Vq 

where q = Pw + ghn, 

Pvv G H^{0,) : -Apvv = trace[(I)u)^] Pvv\g = 

hn€ H^{n) : -Ahn = hn\s = h 

This splitting of q is natural since pw is the Lagrange multiplier of the volume preserving condition 
(divu = 0), and yV/i^^ is the variational derivative of the potential energy = J ~ f gzdxdz. 

nn{z>o} ff^nizKO} 

Proposition 2.1. The conservation of energy for (E-BC) system is 

E= I -\v\'^dxdz + y = I -ipAftp + -gh'^ {I + \Vh\^\yh dS = constant. 

Proof. The fact that E is conserved is a consequence of the Lagrangian formulation of the Euler 
equations. Here we will present a direct proof. Differentiating with respect to t 
dE f dv , , f I. 



[ v-^dxdz+ [ -\v\^v ■ N dS + g [ hv ■ N dS 
Jn dt Js 2 Js 

[ V -V— + v - Vqdxdz + [ -\v\^v ■ N dS + g [ hv ■ N dS 
Jn 2 Js 2 Js 



dt 



qvNdS + g / hvNdS = G. 
s Js 

□ 

Note that the conservation of energy implies bounds on Af'^ip and on h in L'^{S). Thus to 
derive higher energy estimates on ip and h we need to relate the Dirichlet to Neumann operator 
M to differentiation on S. This is done by standard estimate on harmonic functions on O which 
imply that as far as estimates are concerned, N behaves like ^— plus lower order terms. These 
estimates were carried out in jSZ06) without detailing explicitly how the constants in the inequalities 
depend on the various norms of the solution. Here we will detail how the constants in the energy 
inequalities depend on several derivatives of the solution in the L°° norm. For the remainder of 
this section we assume that S is given by the graph of a smooth function h and that \\h\\^r5,p^^2•J < eq 
for some 2 < p < A, and Eq sufficiently small. Thus for any £ > we do not need to distinguish 
between H^{S) and H^{M?), i.e. for any (p G H^{S) we have 

W^Wn'^is) + \\h\\HHs) ~ ll95|b<'(R2) + 
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Proposition 2.2. For any (p defined on S 

(2.1) \\Af^\\w2,^(R2) < \\D^\\w3,c^ + WA^ifWi^ 

The proof of this proposition follows from standard arguments applied to double layer potential 
and will be given in appendix [El Note that inequality (12. ip can be improved by using norms, 
however (12. ip is sufficient for our purpose. One should compare the above proposition to the 
detailed behavior of as an analytic function of h obtained in [CNOO] . 

Based on this proposition and the work done in [S Z06j we make the following observations : 

• Since Ash = {l + \Vh\^y^[{l + \d2h\^)dfh-2dihd2hdl2h+il + \dih\^)d^h] and the mean curvature 
is given by k = —(1 + |V/ip)2 A5/1, then for A; > ^ > 2 

\\h\\le^^2)^ J + hi-A's)hdS ^ J h^ + Ki-AsY-^KdS 

• By defining 

||¥'|Ih«+1/2(5) = ||'^|Ih«(R2) + \\A2'D^ip\\l2(^2) 

then from the single layer potential formula (IF.3h it is easy to see that 

lly'lUaCS) + II^V'lUaCRa) ~ llv^llLSCIRa) + ||AAv3||2,2(ir2) 

and by interpolation 

\\ip\\Hi+i/2(^S) ~ ll¥'ll//^(5) + W-^^T^^vWl^iR^) 

• For any function g defined on 0, and for any function ip defined on S 

(2.2) D^ip'n{N,N) = M^{ip) - 2V N{-A)-^{DNn ■ D^^n) - -^fiN) ■ {Af{^)N + V^^) 

(2.3) Ag = Asg + nVN9 + D^g{N,N) on S ^ 

(2.4) (-As - Af^)^ = kN{<p) - 2Vn{-A)-^{DNu ■ D^^n) - AA(iV) • {Af{^)N + V^^) 

Consequently from equation (|2.4p M'^ is like — plus lower order terms (involving derivatives of 
h), since after integrating by parts twice we have 

[ ^{-Asip)-\Mip\'^ dS = [ K^M{ip) + ^M{N)-{M{ip)N + V~^^) dS-2 [ D Nn{V ^h) dx . 
Js Js Jn 

Therefore for i >2 

11^111^^(5) + II^IIh^ ~ / \^\' + \^f'v>\'ds+\\h\\l,. 

J S 

• Since the Euler flow is assumed to be irrotational, then from [SZOBj (equation (3.6)) 

(2.5) DtK = -{AsV^Pn) ■N-2U- {{D^\sWn), 
and from proposition 4.3 of [SZ06], we have 

Proposition 2.3. ||VV'-hII_h'<'-i/2(5) ^ \\DtK\\fje-5/2(^g-^ + \\h\\fje-i/2 + VE, for £ > 3. 

Sketch of the proof. Since the above statement is contained in the afore cited proposition we will 
only present the idea of the proof: on S, split v = N-Wipji into tangential and normal components, 
compute the surface divergence and curl of these quantities, and use ()2.5p 

-(As VVw) ■N = DtK + 2U- {{D'^ViIjh) 
to obtain the stated estimate. 
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Let and denote the tangential and normal components of v. The tangential divergence of 
i>~^ is given b}0 

(2.6) V-u^ = {AsVtl;H)-N + {D'^V^pH)-'n = -DtK + 0{Y^ \D^^\ + \D^ h\f). 
and the tangential curl by 

(2.7) LoJiXi) • X2 = n(Xi) • {X2 • vvvw) - n(X2) • (Xi • vv^-h) = o{J2 I^'V'I + \D^h\f. 

i=i,2 

where {Xi,X2} an orthonormal frame of TS. Thus 

(2.8) \\Du'^\\jji-5/2 < \\DtK\\Hi-5/2(^s) + + \\Dh\\Hi-3/2 + \/^. 

To bound let = N-^ ■ Wipji and note that 

Mu^ = VNiVNn (VV'w) • Nn) - A-^[A{VNn (Ww) • Nn)] 

V-v'^ = V -v-Kv^ - Vjv(Vjv« (VVw • iVw) + iV • I?(VV^H)(AA(iV)). 

Since V • P = D{Vipn) ■ DNn, we obtain Mi^^ = -V ■ v'^ + 0(Ej=i,2 l-^^V'l + \D^h\f ) which 
together with ()2.8p imply 

IIDA^-VV^/^T^ 11^^-5/2(5.) = \\Dv\\hI-5/2(^s^ < \\DtK\\fjl~5/2(s^ + \\D'4)\\fjl-3/2(s^ + \\Dh\\Ht-3/2(s)+^- 

Since N is equivalent to one derivative in norm, we obtain the stated bound. □ 

These observations imply 
Proposition 2.4. Assume that v = Vip-H and h solve the (E-BC) system, then for £ > 3 

\\Ht)\\m{s) + lb(i)llH*-i/2(5) ~ \\DtK{t)\\fji-5/2(^s) + Mi)\\m-^{s) + 

Commutators estimates. From [SZ06] we have for any function / defined on 0, and ip defined 
on S, and where S is moving by the normal component of the velocity v = Wip-H 

(2.9) AV/ = VA/-(I?'V'w)(V/), 

(2.10) Dt^n = {Dt^)n + A-i(22?Vw • ^Vh) 

(2.11) AA-V = A-^Dtf + A-\2D^^n ■ D^A-'f). 
These equations lead to the following commutators formulas 

(2.12) [As, Dt]f = 2V^f • {{D^lTsWn) + V^/ • AsV%bn - KV^TfV^n • N 

(2.13) [Dt,M]f = VN^-^2D^^n ■ D^fn - 2D^i^n{yfn,N) + Vfn • ND^i^n{N,N). 

which together with (j2.4p imply that for s > 2 these commutators are bounded operators on spaces 
given in equations (4.23) and (A. 14) of [SZ06] . Amplifications of these bounds are given in the 
proposition below. Recall that u = h + iA^^'^tp. 

Proposition 2.5. Assume thatv = V'^'H ctnd h solve the (E-BC) system. Letw{t) = {A^h{t),v{t)) 
defined on S, then the following commutator estimates hold 

(2.14) \\[As,Dt]f{t)\\H^^s) < Mt)\\w^>^\\fmH^+HS) + \\fmw^>^\\w{t)\\Hi+^s) 

(2.15) \\[M,Dt]f{t)\\Hi^s) < hWlli^2-ll/(i)ll//^+i{S) + \\fmw^.oo\\w{t)\\Hi+2^s) 

(2.16) \\[As,M]f{t)\\He^s) ^ lk(t)||iy4,^||/(t)||/^.+2(5) + ||/(t)||ty3,^|k(t)||j,.+5/2(s) 

^Here we introduced the notation A = 0{B) to mean \\A\\h- < \\B\\h- ■ Thus A = 0{Y.^^^ ^ \D^ ip\ + \D^h\f) 
imply \\A\\h^ < (IIDV^IIwi.- + \\Dh\\^^,^){\\DnHs + WDhUA- 
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Proof. From equation (j2.12p we note that [As,Dt] is an operator of order 2 with coefficients 
depending on second derivatives of h and Vipn on S. Thus by Holder and Sobolev inequahties we 
conclude 

\\[As,Dt]fmH^iS) ^ iWVMw^,^ + \\h\\w^oo)\\f{t)\\jj,+2 + \\f{t)\\w2,oo{\\VMH^+2 + \\h\\Hi+2), 

and by proposition 12.21 we conclude inequality (|2.14p . The proof of (|2.15p differs from the above 
by the way we treat w = A~^2L'^V''W " D^fn- This is done by applying vector fields Xa, which are 
tangential to S, to the equation 

Aw = 20"^ ^-H ■ D'^fn w\s = Q 
to obtain ||VX^tD||2,2 bounds. The proof of p.l6p is similar to the proof of p.l5p . 

Remark. Proposition \2. 51 holds with I replaced by £+^. This follows from the definition of H^~^2 . 

Equation for k and energy estimates. In the presence of gravity the evolution of k can be 
derived from |SZ06| and [SZ09] section 6. problem II. 

Proposition 2.6. The evolution equation of the mean curvature k is given by 
(2.17) Dfn + {-VnpWk = R 

where R = 0(^o {D^V^pnl^ + \D^h\^) m norm. 

Proof. From |SZ06| equation (3.15) we have 
(2.18) 

D^^K = -DtAsv • iV - 2n • {D'^\sDtv) + A^VV^w • {D^^h){N)^ + 2[P (((L>Vh)(A^))^) 
+ Ii{{D'^\sVijHV)] ■ {D^\sV^h) + 2n • {VVMs? - 2{{D^nT{N)y ■ Il{{D^nTNy) 
From the equation for [Dt, Ag] (I2.14p and Euler's equations Dtv = —Vp — ge^ we obtain 

DfK, = N ■ AsVp + R 
where R = 0(J2l iD^V^pnl"^ + |L>J7i|2) in norm. Computing N ■ A^Vp 
N ■ AsVp =N ■ AVp + N ■ {AsVp - AVp) 

=V N^P - N ■ (kV N^P + D^{Vp){N, N) 
=VN^P-Nn • {KnVNn^P + D\Vp){Nn,Nn)) 
=Vn^P - ^Ninn^NnP + D'^p{Nn,Nn)) 

+ VnpMk + kVp • VnN-h + 2D'^p{N,VNNn)- 
Since Ap = — tr(L'v)^ = — iAjV^/'T^p and p\s = 0, then 

\\l^H'^NnP + D^p{Nn,Nn)\\He{S) = W^P- '^SP\\h'^(s) = \\^^{Dvf\\HHs) 

l|A(KwVAr^P + D'^p{Nn, Nn))\\Hi~3/2(n) < W^^iDuff^t+^^^f^^^ Mm-^s) 

which implies that N ■ AsVp = VnpJ^k + R where R = 0{^l {D^V^Jnl^ + l^-'^P)- This gives 
the stated evolution equation for k. □ 

Based on this equation we define the high energy Ek(t) = ekdS + E for > 3 as 

, , [ ekdS= [[M{-As)''DtK]{-As)''DtK+{-VNP)m-As)''K\US 

(2.19) Js Js 

={Mi-As)''DtK, {-AsfDt^^) + {i-VNP)M{-As)''K,M{-As)'K) 
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where (,) denotes the inner product on L^(5). Note that since 

q = Pv,v + ghn -Vnp = -^nq + gN ■ 63 > g - ceq 

then Ek ~ ||V^7^ || ^3^+5/3 + | |/i|| ^2^+3(5) • 
Proposition 2.7. For t > 2, 

(2.20) Efc(t) < Ek{2) + ^ ||u(s)||H/4,.o^fc(s) ds 

Proof. To compute -DfCfc we proceed as follows: 

• \DtV iyp\L°° < |M(t)|iy4,oo. This follows from the identity Vp = Vpv,v + V/it^ — 563 and the fact 
that 

AA-V = A-iA/ + A-i(2I)Vw^'A-V). 

• AAA(A5)*^K = AAA(-A)'=K+ [Z?i,AA](-A5)'=K. From (IXT?]1 we have 

|([A,AA](-A5)'=ac,AA(-A5)'^k)| < \u\yv^.o. Ek{t). 

k-l 

. AAA(-A5)'=K = J^{-As)''DtK + ^ AA(-A5)nA, -A5](-A5)'=-^-^At. Since [Dt, As] is a sec- 

i=0 

ond order operator we have from (j2.14p 

{AfDt{-As)''K,Af{-As)''K) = {M{-AsfDtK,M{-AsfK)+0{\u\w,,^Ek{t)). 

• Computing dt{M{-As)''DtK, {-As)^DtK) 

dt{M{-As)'DtK,i-As)''DtK) = 2{Dt{-AsfDttiM{-AsfDtK) 

+ {[DtM]{-AsfDtK, {-AsfDtK) 

Since [Dt^N] is a first order operator we have from (j2.15p 

\{[Dt,M]{-AsfDtK,{-AsfDtK)\ < \u\wi.^Ek{t) 
fc-i 

. Dt{- AsfDtK = {-AsfD'in + J^(-As)'[A, -As\{-Asf-'-^ DtK. Since [Dt,N] is a second 

^^^^ i=0 

order operator we have by (j2.14|) 

{Dt{-AsfDtKM{-AsfDtK) = {{-Asf D'iKM {-AsfDtK) +0{\uy^,^Ek{t)) 
Using equation (j2.17p we obtain 

j^Ek = 2{{-Vp)N{-AsfK - {-Asf{VNP)MK , M{- AsfDtK) + O{\u\w,,o. Ek{t)) 

Commuting Vj^p and J\f with {—Asf and using the fact that [AA, A^] is a second order operator 
with error bounds given in (|2.16p we conclude that 

^Ek{t) = Oi\u\w'i,o.Ekit)) 

and thus 

Ekit)<Eki2) + j C{eo)\u\w^,^Ek{t)ds. 

□ 
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Thus from the assumption ||/i||^y5,p(]{52) < Eq, reference |SZ09j . and a limiting argument to go 
from smooth initial data to data in H'^ where = 2A; + 3 ^ 5/2, we have local well posedness 
and the energy estimate stated in the above proposition. 

3. Space time resonances of quadratic and cubic terms 

Expanding (l^/'VF) in powers of h and ip, setting g = 1, and keeping track of quadratic and cubic 
terms we obtain 

( dth = Aij-V- (hVi;) - A{hAi;) - \ {A{h^A^^\)) + A^{h^A'i\)) - 2A(hA{hA^))) + i?i 
^•^■^^ \ dti) = -h- i|V^|2 + \\AiP\'^ + AiPihA^ij - A{hAiP)) + i?2. 

where Ri and R2 are of order 4. We refer to the book of Sulem and Sulem |SS99j for the above 
expansion (also see the remark at the end of appendix [F]) . 

Writing the equation in Fourier space. Recall that 

u"^ {h + iA\^l:) , uo=^e-'^''''\ho + iAho) and f e''^''\ = e^'^''\h + zA'/^^). 
Writing Duhamel formula for / in Fourier space yields 

(3.2) f{t,0 = uo{0+ E /7e'^'^^^"^"^.(^,^)/^i(^,^)/^2(^>e-r/)d7?ds 

Tl,2=±J = l ^ 



+ E E^J'^1'^2'^3 / jje'"''^^'^'^'^3mj{^,r],a)f-ri{s,v)f-T2{s,(T)f-r3is,(-ri-a)dr]dads 

Tl,2,3=± i=3 ^ 

rt 



where Cj^±^± and Cj^±^±^± are complex coefficients, /+ = f, f- = f, and R = Ri + iA'^R2 is the 
remainder term of order 4. The phases are given by 

'/'±,±(e,r/) = |Cl'/'±|r/|^/'±|e-r/|i/2 
^±,±Ai^V,^) = lei'/' ± H^'^ ± kl'/' ± IC - r? - ct\^'\ 
and the multilinear symbols are defined by 

def 1 

= 2 |^|l/2|^_^|l/2 + 

m3(e, r/, a) (le - - ct|3/2 + i^ii^ _ ^ _ ^|i/2 _ 2|^ - - - a\^'^ 

-4(e, r/, a) [\i-ri-a\'/'-\i-v\\i-V- ^1^/^) • 

Note that mi and m2 are homogeneous of degree 3/2 and that and 771,4 are homogeneous 
of degree 5/2. Also, note that since these multilinear forms are homogeneous, we only need to 
estimate them on the sphere + |r/p = 1, or + \r]\'^ + |crp = 1 and extend all estimates by 
homogeneity. Moreover the exact form of the above equation is not really important; thus in order 
to focus on the information which is relevant to us, we shall ignore from now on the distinction 
between /+ and whenever this notation occurs. 
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Examination of the quadratic symbols. The symbols rrii and m2 have two important features, 
they vanish when one of the Fourier coordinates rj oi — rj) is zero, and they are not smooth. 
These two facts are made more precise in the following lines. 

Notice that the vanishing of mi and m2 is crucial: as we will see, it corresponds to a null property 
on the time resonant set; on the other hand, the lack of smoothness is a hindrance, since it prevents 
one from applying the standard Coifman- Meyer theorem [CM78] . 

We always use the convention that A stands for a smooth function in all its arguments and start 
with mi: 

. If |??| « 1^1 ~ 1, mi(e,r?) = \^\y^A{\v\'/^^,^). 

. If|C|«|7?|~l, mi(C,r?) = |e|^(|,r7). 

. If - v\ « 1^1 ~ 1, V) = \^- (l? - e) . 

Now m2- 

. If \v\ « m ~ 1, m2{(,v) = \r]\'/'A (|^|i/2,i,^). 
. If|C|«k|~l, m2{^,v) = \e^^A{v,0- 

. If - r?| « 1^1 ~ 1, m2(e, r/) = le - vl'^^A (|e - v\'^^ ^,^) ■ 

Space and time resonances. We define the bilinear and trilinear time resonant sets as 

^±,± = {0±,± = 0}, ^±,±,± = {cP±,±,± = 0} 

respectively. We define the bilinear and trilinear space resonant sets as 

^±,± = {v^</.±,± = 0}, ^±,±,± = {V^,,(/.±,±,± = 0} 

respectively. The bilinear and trilinear space-time resonant sets are given by 

^±,± = y±,±n^±,± , ^±,±,± = ^±,±,± n 5±,±,±. 

respectively. 

Examination of the quadratic phases. The phase = |^|^''^ + Ivl^^'^ + ~ is better 
behaved than the others, since it only vanishes at {(,,1]) = (0,0). Therefore we focus on the three 

other quadratic phases, namely (j) j 4' h 4'-\ Up to multiplication by —1, and permutation 

of the three Fourier variables rj, — r], these three phases are the same. Let us consider 

A small computation shows that 

£r__ = {rj = or ^ - r/ = 0}, 
and the vanishing of (p may be described as follows 

(3.3) if « ~ 1, <A-(e,r/) = |?7l'/'^(hl'/',^,e) with ^(0, •, •) = -1 

(the case ^ — rj = being identical up to an obvious change of variables). Finally, along the surface 
= 0}, (p does not vanish, but is not smooth: 



(3.4) if lei « \v\ ~ 1, 0-(e,^) = A (|r/|V2, + |^|i/2^' (^|^|i/2, 

with ^(0, •,•) = -2|7/|i/2. As we shall see, it turns out that quadratic terms can be treated simply 
by a normal form transform, thus there is no need to investigate the quadratic space resonant set. 
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Examination of the cubic phases. Cubic phases fall into two categories: some have relatively 
few time resonances, we call them 'weakly resonant phases'; and some give a large space-time 
resonant set, we call them 'strongly resonant phases'. 

The weakly resonant phases: (p h+j ^Pa h> ; 4* As in the case of quadratic phases, 

the phase 0+++ = I'^l^''^ + Ivl^^"^ + l^l^^"^ + 1^ ~ ~ ct|^/^ is easily dealt with, since it only vanishes 
at iC,rj,a) = (0,0,0). 

The four other quadratic phases, cp h+) 'Pa hi 0+H > ! identical, up to multiplication 

by —1, and permutation of the four Fourier variables rj, a, ^ — rj — a. Let us therefore focus on 

It is easily seen that 

= {rj = a = or a = ^ — rj — a = or r] = ^ — rj — a = 0}. 

The vanishing of on this set can be more precisely described as follows 

(3.5) if |r/|,|a| « ~1, <P = -|7?|i/2-|a|^/2^^(e,7?,cT)[r? + cT]. 

(recall A{^, r], a) [r] + a] stands for a function smooth in its three first arguments and linear in the 
fourth) . The other cases a = ^, — r] — a = and t] = ^ — r] — a = are the same up to a change of 

variables. Finally, (j) is not smooth along the axes = 0} |J{?? = 0} \J{(7 = 0} \J{C — r] — a = 0}. 

In a neighborhood of these axes it can be written in a form similar to (j3.4p . 

As for the quadratic phases, it turns out that a normal form transform is sufficient to treat the 
weakly resonant cubic terms, therefore, we do not investigate the space resonant set. 

The strongly resonant phases: <j) \., (p \ , Once again, these three phases are identical 

up to a permutation of the Fourier variables r], a and ^ — rj — a. We therefore focus on one of them, 
namely 

'i>-+ = \i\"'-H"'-W\"' + \i-v-a\'". 

In this case, the time resonant ^z'. |_ set has dimension 5, and (j) ^ vanishes at order 1 on it: this 

makes a normal form transform almost necessarily unbounded. We therefore need to turn to the 
space resonant set, which is 

^-+ = {e = r? = ^}, 

and this set is contained in 3/^. |_. In other words |_ = M y and it has dimension equal to 2 

and the phase vanishes at order 1, therefore an argument based on an integration by parts in the 
Tj and a seems doomed to fail. 

The way out of this problem appears if one develops d^cj) ^ and daCp \- in a neigh- 
borhood of 5^ +. Namely, if \r] - i\,\(T - i\ « |^| ~ 1, 

dr,4-+ = ((^ " + " ■ il il + ''^^^'^ - ^ - 0, - e, ^ - 0] 

d,<P..+ = (^(^ _ ^) + ^(^ _ ,,) . + A{i, V, <t) [(r? - e, ^ - 6, (^ - - 0]) 

55</._„+ = ((2^ - ^ - ^) + ^(^ + ^ - 26 • il il + ^, ^) [{r]-t'J-0,iv-C,'T- O) 

(in the above expressions, A is smooth in the three first arguments, and bilinear in the arguments 
between brackets). The above expressions imply that d^cj) = —dr}(p — da(j) up to second order terms. 
Therefore, if — CI) I''" ~ CI << ICI ~ 1) 

(3.6) d£_4>—+ = A{i,r],a)[dj^4>^_+,d„<j)—+\. 

This identity will enable us to convert C derivatives of (f) into r] and a derivatives of (f). The former 
occur when applying to our multilinear expression (which corresponds to the x weight), and are 
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problematic since they come with a t or s factor; the latter are harmless: an integration by parts 
gets rid of them. See Section [7] for the details. 



4. Normal form transform 
Integrate by parts in s the quadratic terms of (13. 2|) with the help of the formula ^ dgc'^^'^'^''^ = 
g«s0±,±_ Doing so, the dg derivative will hit f{s,ri) or f{s,^ — rj); then use ()3.2p to substitute for 
ds f{s,ri) or ds f{s,(, - rj). This gives 

e»^'^n,-2 mii^, 7])f{s, r?)/(s, ^ - rj) drids 



f{s,r]) f{s,(, - r])dr] 



9. 



f{s,v) f{s,(.- V) 



drjds 



2 fi,r2=±i = l,2 



) fis, V) fis, cr) fis, ^-'n-cr)dr]dads + 




ism'^^+riM'^-') ri) j^^^ ^^Q^^ _ ^ 1 +{ symmetric terms }, 

where the "symmetric terms" come from the fact that ds may hit either f(s, ^ — rj) or /(s, rj), thus 
these symmetric terms look very much like the above ones with rj and — rj exchanged. The term 
Q corresponds to terms of order 3 and higher in e~^^^^^^ dtf i.e. 

4 

def 



l,T\,T2,TZ 

Tl,2,3=± J=3 



"^jiii o')u{t, r])u{t, a)u{t, ^ — i] — a) dr] da + R{t,(,). 



The normal form transform that we just performed gives quadratic (without time integration), 
cubic, and higher order terms. Also, cubic and higher order terms occurring in (|3.2p need to be 
taken into account. In the end, we see that / can be written as a sum of quadratic terms of the 
type 



(4.1) 



/ 



e^*^±'±^^^/(t,r?)/(t,e-^)dr? 



with / = 1,2 



(here it is understood that the symbol it stands each time for either + or — ), cubic terms of the 
type 




H^lj^ _ ^)/(^^ ^)/(^^ ^)/(^^ ^-r,-a)drjdads with /, j = 1, 2 



(4.2a) 

(4.2b) J^Jj e'"''^-^'^mi{^,ri,a)f{s,rj)f{s,a)f{s,^-r]-a)d'ndads with/ = 3,4, 

and terms of order 4 of the type 



(4.3a) 
(4.3b) 



1/2 



R{s,Ods 



^(l^l^^^+-l''l^^^)^^/(.,.)Q(.,e-.)d.d. 
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5. Estimates for the quadratic terms 

In this section, we derive the desired estimates in L°° and L^{x^dx) for terms of the form (j4.ip . 
We shall denote a generic term of this type by 

(5.1) 9i{t,0 '= I e'"^f,{^,v)f{t,v)ht,^-v)drj, 

where 

9 = (P±,± and n{^,r]) = — , 

with the index / equal to 1 or 2, and it either + or — . 

Of course, we adopt this lighter notation because the precise value of it or Z will not affect the 
argument which follows. 

Preliminary observations and reductions. Let us first consider the symbol rj) occurring 
in the definition of gi: regardless of the precise indices / and ±, it, it is 

• Homogeneous of degree 1 . 

• Smooth if none of ??, ^, ^ — vanish. 

. Of the form ^(C, r?) = ^ (|7?| V2, ^, ^) if |^| « |^| ~ 1. 

. Oftheform;u(^,7?) =^(|^-r?|i/2,|52.,^) if 1^ - ??| « ICI ~ 1- 
. Of the form ;u(^,7?) = \C\'/^A ||,??) if iCl « \v\ ~ 1, 

the last three points follow from the developments given in Section [3l Thus we conclude that /i 
belongs to the class Bi, see appendix [C] for the definition. 



Next, define a cut off function xiC^v) which is valued in [0, 1], homogeneous of degree 0, smooth 
outside of (0,0), and such that x(C)^) = in a neighborhood of {r/ = 0}, and xiC^v) = 1 in a 
neighborhood of — ?? = 0} on the sphere. Then one can split gi as follows 

+ J [I - X{C,V)] f^iC,r])f{t,r,)f{t,^ - ri)dv. 

By symmetry it suffices to consider the first term of the above right-hand side, which corresponds 
to a region where \r]\ > |^|, |^ — 7?| (the interest of that condition is that ^ derivatives always hit 
/(^ — rj), which corresponds to low frequencies). Thus in the following, we shall consider that 

9i{0 = I e''^x{^,v)KC,v)f{t,v)f{t,C-v)dv 
which means in physical space, using the notation introduced in appendix [Bj 

g, = e'''''''B^,{n,u). 

Actually, depending on the it, it appearing in (/>±^± the u can be u or u; but as indicated in Sectional 
we ignore for the sake of simplicity in the notations this distinction. Also, it is important to notice 
here that, since fi £ Bi and x ^ t^X £ -^i- 
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Quantities controlled by the X norm. Before estimating the term gi in X, we give some general 
estimates that will be useful in the whole proof. Interpolating between the different components of 
the X norm gives the following lemma. Let Pj, P>j, P<j denote the Littlewood-Paley projections 
defined in Section [Al 



Lemma 5.1. (i) If 2 < p < oo and k < N + ^ - 1, 



VPyjU 



< 2 



t \Mx- 



(ii)If2<p<oo andk < N +^-1, \\V''u\\p < t 

(ttt)Ifl<p<2, ||/||p<i(^"')'||n|U. 
Proof. The proof of (i) is a standard application of (jA.ip and (|A.2p : 



mix- 



VP>jU 



^(-V+fc+l--) II II 



^ oj(-^+fc+l--)ll II ^ oi(-^+fc+l-^),5|| II 



For (iz), we estimate separately low and high frequencies: 

II V^ulL < \\P<jV^u\U + ll-P>,-V'=n|L < 2^''=||n|L + 2- 



<2^-^t-^+l||n||x + 2- 



>] V "-lip rO ^ II "-IIP 
2 



j(_7V+fc+l_|) X 



Optimizing the above inequality over j gives the desired conclusion. Finally, {Hi) follows from 
interpolating LP spaces between weighted spaces. □ 

Bound for V'^e~**^^''^5i in with < A; < 4. The idea for this estimate, as for many estimates 
which follow, is to use Sobolev embedding in order to make sure that Lebesgue indices are finite 
when applying Theorem lC.il 

Using successively Sobolev embedding. Theorem IC . 1 1 and Lemma ISTT] one gets 



< 



< IM'^+i 



^llvyi'16 IfIIi6 



< 



M'x- 



Bound for xgi in L^. In Fourier space, xgi reads 
(5.2a) r{xgi) =8^ 

(5.2b) 

(5.2c) 

(5.2d) 



e''^x{i.r,)i,{7^,i)f{t,il)f{t.i-il)dr) 

+ j e'"^d^[xiC,v)l^{r],0]fit,v)f{t,C - V)dr] 
+ [ e''hd^4>{^,r,)x{^,v)f^{ri,0fit,v)fit,^ - V)dv- 



Bound for / f 5. 2b]) in L? . Using successively Sobolev embedding. Theorem IC.ll and Lemma I5.lt one 
gets 



B.^,[u,e^''^''\f 



< ||yln||,^i,4 ||x/||2 < ||n| 



< 

2 ~ 
2 

X- 



B^,[u,e^''^''\f 



1 4 
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Bound for 115. 2c]) in L^. Let us take a closer look at xCCi 0- This is a symbol in Bi, which 

furthermore vanishes at order ^ in ^ = 0. Therefore, by Lemma |C.3[ 

[xiC, V)K^, 0] = f^i + j^|17^/^2 + 1^ _ ^| At3 with (/ii, /i2, /^s) G So X ^1/2 x 

Thus 



+ 



+ 



1 



1 



(5.3a) 
(5.3b) 

(5.30 . ,^ 

The term ()5.3ap is easily estimated, so we skip it, and consider next the term ()5.3bp 

B^^{u,u) < ||S^2(n,n)||4/3 < 



r^^3iv,0fit,v)fit,^-v)dr]. 



KlObDii, 



^1/2 



< 



\U\ 



Finally, we need to estimate the term (j5.3cp : by Theorem IC.ll and Lemma lA.ll 

11/2 1 



< 



Uu\\Jf\\,/,<\\u\ 



< 



|2 

\x- 



Bound for 15. 2d\) in L^. Since 



" 1 |g|3/2 ^ I |g->7p/^ ' keeping in mind that ji vanishes at 
order 1/2 in ^, the symbol appearing in ()5.2dp can be written 

def 1 ~ ~ 

d^Ht ??)x(C, ^)^(??, 6 = fJ-iit V) + 1^ _^|i/2 ^2(^, rf) with (;Ui, ^2) G -61/2 x fii. 

We simply show how to estimate the term associated to the symbol |^_.^|i/2 ^2(C; v)-, the symbol 
A*i(?)'7) being easier to treat. Using successively Theorem IC.ll Lemma lA.ll and Lemma l5.ll one 
gets 



|^-ry|i/2 



fJ'2 f{t,v) fit,C- V)dr] 



<t\\M\^ 



<t\\Au\\. 



yll/2- 



1-2^0 



-l+25o + iv+23^[l-2<5o+5]^(i_25o)5 



< tt 



<A\u\\l. 



\x 



where > denotes a small constant. Notice that the last inequality holds since N and 5q have 
been picked to be large and small enough respectively . 

6. Estimates for the weakly resonant cubic terms 

In this section, we derive the desired estimates in L°° and Lp'{x^dx) for terms of the form 
either (j4.2ap or ()4.2bp corresponding to weakly resonant phases as defined in Section [3l We shall 
denote a generic term of this type by 

ft 



(6.1) 



52(e,t) 



def 



e*"^M(C, r], a)f{s, a)f{s, r])f{s, ^ - r] - a) dr] da ds, 
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where 

4> = 4>+++ or or (/)+_+ or or (j) 

and 

l^{£.,V,(^) =m{C,V,(^) or — -——mj{^-r],a) or — — -mj{^ - ri,a), 

with the indices i equal to 3 or 4, j, k equal to 1 or 2, and it equal to + or — . 

Of course, we adopt this lighter notation because the precise form of (p or fj, will not affect the 
argument which follows. 

Preliminary observations and reductions. We begin with the symbol //(^, rj, a) occurring in 
the definition of §2 ■ regardless of its precise form. From Section [3] and the observations at the 
beginning of Section [5] we have 

• It is homogeneous of degree 5/2 in (^,77, a). 

• It is smooth except on = 0} U {/? = 0} U {a = 0} U - ry - cj = 0} U - r/ = 0}. 

• It might have a singularity for ^ — 77 = of type, in the worst possible case, |^ — ryl^/^. 

• It vanishes to order (at least) 1/2 in ^. 

• It belongs to the class 7^/2 (see appendix [D] for the definition). 

Next, define cut-off functions xi; X2, Xs such that 

• Xi, X2, X3 are valued in [0, 1] and xi + X2 + X3 = 1- 

• Xi) X2, X3 are homogeneous of degree and smooth outside of (0, 0, 0). 

• Xii^i 7?, t) = in a neighborhood of {(7 = 0}. 

• X2iS,, rj, a) = in a. neighborhood of {r/ = 0}. 

• X3(Ci rj, a) = in a neighborhood of — r/ — cj = 0}. 
Then one can split g2 as follows 

^2(^) = YlJ^JJ e'''''Xk{^,V,(^)l^{^,ri,(^)f{s,^)fis,r])f{s,C-r]-a)dr]dads. 

By symmetry it suffices to consider the first summand {k = 1) of the above right-hand side, which 
corresponds to the region \a\ > |^|, |r/|, \S,—ri — a\ (the interest of that condition is that ^ derivatives 
always hit /(^ — rj — a), which corresponds thus to low frequencies). Thus in the following, we shall 
consider that 

92(0 = J^JJ ^'"^XiiC, V, f^)/^(C, V, (^)f{s, v)f{s, (T)f{s, ^-r]-a)dT]da ds. 
Since xi ^ '^0 and fi £ '^5/21 we will constantly use in the estimates which follow that xi^ ^ '^/2- 

Finally, we will need the following lemma 
Lemma 6.1. If (f> is either (f>^-^-^, (j) |_-|-, 0-| |_, , or 4> , then ^ can he written as 

^ = "° + W^''' + W^"' + + 

where (z/q, 1^1,1^2, 1^3, i^i) e '^'-1/2 x Tq x Tq x 7^ x T^j. 

Proof. Relies on (jS.Sp . □ 

Remark. Depending on which one o/(/)++-|-, (j) |-+, (j)^ y, (f)j^^ , or (j) is considered, one of 

the four symbols (z^i, V2^i^3^ ^^4) can he taken equal to 0. 
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Control of e dtf. The normal form transformation, i.e., integrating by part with respect 
to time, introduces terms of the type e~**^^ ^ dtf ■ This explains the importance of the following 
lemma. 

Lemma 6.2. If 2 <p < co and < A; < 10, 



-2+^ + 



Proof. Differentiating ()3.2p with respect to t and applying e * gives that e-^*^'^ 9J IS a sum 
of terms of the type 

(6.3a) j mj{^,r])u{t,r])u{t,(, - r])dr] j 



1,2 



(6.3b) 



"^j ('^1 o')u{t, r])u{t, a)u{t, — rj — a) dr] da j = 3, 4. 



plus remainder terms. With the usual justifications. 



V^(|6.3aj) 


< 

p 




V''(j6.3b[) 


< 

p 





-2H 



2p 



3p 



-3 , , 



\u\ 



Xi 



and the estimate follows by using Proposition IF. II for the remainder terms. 



□ 



Bound for V'^e **^^'^^g2 in , for < < 4. Integrating by parts in s in ()6.ip using the relation 



up * 



52 (t, 6 



e''"^^Xi(C, ^' ^)/^(^' ^' cr)/(s> i - ri - a) drj da 



1 t 




e**'^-7Xi(C, ^, cT)/i(e, ct)5, /(s, a)f{s, v)fis, i - r] - a) 

Kp L 



dr] da ds. 



The boundary term at s = 2 is determined by the initial data and is easy to estimate. Next, replace 
by its decomposition on the right-hand side of ()6.2p . The term vq is harmless of course, and so 

is the term j^^i'i, due to the vanishing of /i in at order 1/2. The three remaining terms can 
be treated in essentially the same way, thus for the sake of illustration we only retain the term 
l^|i/2 ^2{^, Vi (^)- Furthermore, if one distributes the dg derivative occurring in the second summand 
of the above right-hand side, it suffices to consider one of the three resulting terms; we choose to 



consider the case where dg hits f{s,a). Thus J^V e 
the type 

(6.4a) 



k„-itA^/^ 



§2 can be written as a sum of terms of 



da 



(6.4b) 




C'' j^^'^2XiA^S(t, a)u{t, r])u{t, ^-r]-a)dr] 

^A:g-.{t-.)|C| V2 _l_^^^^^^±.slH J(^^ a)u(s, r?)n(s, i-r^-a)dr^da ds. 



Using successively Sobolev embedding, Theorem ID. II and Lemma lA. 11 we get 



1 



n, u 



< 



V B^2tJ.xi [ ^' ^1/2 



n, u 



< 



.11/2 



n 



1 24 ^ 



24 



y^5/2 + A 



W 



,24 IFII24/7 IPI 



24 



\u\\x- 
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In order to estimate ()6.4bp . we split it into two pieces + J^_^■ The first summand is estimated 
using the dispersive estimate of Lemma lA. II and point (ii) of Theorem lD.il 



J" 



-1 



< 



< 



f-i 



ds 



t - s 

*-i ds 



t - s 



t - s 



\u\\2 II^IU^S < -\\ufx. 



The second summand is estimated via Sobolev embedding: 

"* 1/2 1 



J" 



-1 



t 



__^2Xi/ie±*'l'^l'''5 J(s, a)u{s, v)u{s, ^ - r, - a) drj da ds 



\V\ 



(6.5) 



< 



< 



< 



t-i 
t 



t-i 
t 



t-i 



ds 



1 



.11/2 



u 



l^lls ^'^ 



dsf 



WuWal-x \\u\\o ds < -\\u\\\. 
H/2>4 " "8/3 II lis ~ t II 11-^ 



Bound for xg2 in L?. In Fourier space, xg2 reads 



(6.6a) 
(6.6b) 
(6.6c) 
(6.6d) 



+ 



+ 



e''"*'xifJ'f{s, cr)f{s, 'n)f{s, ^ - r] - a) dr] da ds 
'"^Xif^Rs, (J)fis, v)d^fis, i-r]-a)d'qdads 
e*'^a5(xi/i)7(s, a)J{s, v)fis, i - - a) di^ da ds 
e""l's{d^(t))Xif^f{s, a)f{s, r/)/(s, - r] - a) dr] da ds. 



Bound for W. 6h\) in L^. Using successively Sobolev embedding, Theorem ID . 1 1 and Lemma |5.1| one 
gets 

ll(l6lbDll2 = 







) ds 






2 J2 



H/1,4/3 



ds 



< 



Ixfhds < \\ufx- 



Bound for i6.6c\) in L^. Let us look more closely at the symbol xiM- It belongs to T^/^, but it is 
a bit smoother than general symbols of this class: namely, it has (in the worst case) a singularity 
of type — r/l^/^ at ^ — r/ = 0, it vanishes at order 1/2 in ^, and it is smooth for — o" = and 
r] + a = 0. Combining these observations with Lemma ID. 31 we deduce that d^{xifJ-) can be written 



0^ [(Xi^)(^, V, ct)] = / + + ^— + 



|^_^|i/2- l^-r^-a 
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with {ni,iJ,2,H3,fJ'4:) G %/2 X T2 X T2 X T5/2. Thus 



(16:6^ 



V/(5, c^ms, V)f{s, i-'n-a)drjdads 



+ 



+ 




1 



(6.7a) 
(6.7b) 
(6.7c) 

(6.7d) ' I II - I, I 

The term (j6.7ap is easily estimated, thus we skip it, and focus first on ()6.7bp and (j6.7cp . These two 
terms can be estimated in essentially the same way, except that fractional integration is handled 
for the former with Lemma lA.ll and for the latter with Proposition ID.4[ We simply show how to 
deal with (j6.7cp : using Proposition [DA\ 



1 



Yj2f^ f{s, cr)f{s, ri)f{s, ^ - rj - a) di] da ds 
7fi'^f{s, cr)/(s, r/)/(s, i-'q-a)dr]da ds. 



||(I62ED||2< / B 1 ,iu,u,u) 
Finally, the term ()6.7dp is estimated in a very similar way to (|5.3cp . thus we skip it. 

Bound for W. 6d\) in . First compute 
(6.8a) (ICTll -- 



ds <, I ||yl -u-ll^ ||u||4 ||ti||4 < ||n 



3 




1 ^ ^ ^ ^ 
e'^^^s-— t7ttXi/^/(s, cr) f{s, v) f{s, i - rj - a) dr] da ds 



2 |e|3/2' 



(6.8b) 




e''^s\^ — - — ^Xi/"/(s, 'y)I{s, v)fis, i - 1] - dr] da ds. 



2|e-r/-f7|3/2' 



Observe that in the term (j6.8ap the singularity ^ is cancelled by the vanishing of the symbol 

Xi/i in ^. In other words, the symbol ^ |g|3/2 Xi/^ belongs to T2. With the help of Theorem I C.H this 
makes the estimate of ()6.8ap straightforward: 

ft 



|(IM)||2< 



Bi {u,u,u) 

^|3/2 A-lf* 



ds < I s ||yl n||g ||u||g ||n||g ds 



< Il7.lli 



t 2,2 



ss 



s is 3 < t'^llulll 



X 



since has been taken big enough. 

The estimate of (I6.8bp is a little more technical, and uses fractional integration (Lemma lA.ip : 
5q standing for a small constant, 



|(l6:ibDll2<^^ 



< 



1 



\\u\\ 1 



ds 



ds < 



\u\\ 1 



ds 



< llnlli 



l+25o + iv^^(<5-25o+l) l+25o/(|-45o) < ||^||3 .5 



SS 



rv^ II 1 1 

since N has been chosen big enough, and 5o small enough. 
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7. Estimates for the strongly resonant cubic terms 

In this section, we derive the desired estimates in L°° and Lp'{x^dx) for terms of the form 
either (j4.2ap or ()4.2bp corresponding to strongly resonant phases as defined in Section [3l The three 
strongly resonant phases are identical up to a permutation of the Fourier variables. Thus, making 

a change of variables in the r],a integral if needed, we shall only consider the phase (j) h- The 

generic term we consider reads 

(7.1) t) =^ fj j e*^V(C, cj)f{s, rj)f{s, C-rj-a)drjda ds, 

where (p = 4' h ^'I'ld 

H{^,r]) =mi{(„rj,a) or — ——mj{^-ri,a) or — — -mj{^-ri,a), 

with i equal to 3 or 4, j, k equal to 1 or 2, and it equal to + or — . 

Preliminary observations and reductions. Proceeding as in Section [6l and keeping the same 
notations, matters reduce to 

93{t, = fj f 77, ct)/(s, r/)/(5, ^-rj-a)d7ida ds, 

where xi^ e T5/2. 

Bound for xg^ in L^. In Fourier space, xg^ reads 

(7.2a) T{xgs) = j'j j e''^Xl^^f{s, a)f{s, v)dj{s, i - r, - a) dr^ da ds 

(7.2b) + j'j I e'''f'd^{xw)f{s, a)f{s, r?)/(5, ^ - r/ - <t) drj da ds 

(7.2c) + j'j I e''h{d^(l))xMs, a)f{s, v)f{s, ^ - v - <t) dr] da ds. 

The terms (j7.2ap . ()7.2bp and (j7.2cp can be estimated exactly as ()6.6bp . (j6.6cp and ()6.6dp in Section[6l 
This gives, respectively, bounds of 0(1), 0{1), and 0{t^) as t —>■ 00, which yields the desired a 
priori estimate for the norm of gs in L'^{x^dx). To derive the L°° decay of e~^^^^'^ g^ however, it 
will be convenient to prove first that the weighted norm of (73 in L^(x^dx) is bounded as t goes to 
infinity. 

Claim 7.1. The following bound holds: \\{7.2c)\\2 ^ a consequence, \\xg^\\2 ^ H^^Hv- 

In order to estimate (|7.2cp without the time growth factor, the first step is to distinguish between 
high and low frequencies. Denoting 5q for a small constant (whose precise value will be fixed later 
on), we split (j7.2cp as follows: 

(7.3a) = j'j j e^^*5(9g</.)xi/i [l - 6 (^)] f{s, a)f{s, v)f{s, i - - a) dr, da ds 

(7.3b) + j'j j e''h{d^<l^)xilJiQ (^) /(s, a)f{s, v)f{s, ^-r,-a)dvda ds, 

where is defined in Section \^ The point of such a decomposition is of course that there 
only appears in ()7.3bp frequencies (essentially) smaller than s^°, since Supp (^Xi('^i "^i cr)© (?")) 

{n\rilW\<s''}- 
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Bound for (TZSaU. Due to the control in H , it will be easy to estimate ()7.3ap . First recall that 



2|?P/2 + 2|g-r?-cr|3/2 
2|C-r?-(T|3/2 

^ - 7? - 0- 



The singularity in ^ is canceled by /i, thus we can forget about it, replace 



by oi/ Ji'3/2 ™ (|7.3ap , and estimate with the help of Theorem ID.ll and Lemma 15.11 




< 



< 







r 


s 






r 


s 



|e-ry-a|3/2 

(^>5olog(s) 



ds< 



P 



ds < 



In order for the last inequality to hold, we need 



(7.4) 



1 + (5o -A^ + 



13 



< -1. 



Bound for ( 7.3b ). In order to estimate ()7.3bp . a further partition of frequency space is needed. 
Define 

Aj {d,<i, = 0} = {e = ^} 

(7.5) An =^ {9,0 = 0} = {e = a} 

Am =^ {d^cl) = 0} n {d^cl) = 0} = {e = r? = fx}. 
The associated cut-off functions are xij Xii-, Xiih which are taken such that 

• XI, XII, XIII are valued in [0, 1]. 

• XI, Xii, Xiii ai'e homogeneous of degree and smooth outside of (0,0,0). 

• On the sphere \{(,,ri,a)\ = 1, xiii = 1 within a distance of Ajjj, and xiii = if the 
distance to Ajjj is more than Thus on Supp xiii, equation (j3.6p holds. 

• On the sphere |^, r/, o")| = 1, xi (respectively xii) is 1 on a neighborhood of Aj (respectively 
An). 

Then decompose ()7.3bp with the help of these cut-off functions: 
ft 



(7.6a) 
(7.6b) 
(7.6c) 
(7.6d) 





+ 



+ 




e'"Ps{d^(j))xitiQ j /(s, ct)/(s, rj)f{s, ^ - r] - a) d-q da ds 
e'"''s{d^(p)xiXif^Q (^^) f{s, a)f{s, rj)f{s, i - r] - a) dr] da ds 
*^*s(a^(/))x//XiAt0 (^) I{s, cr)f{s, v)f{s, - rj - a) dq da ds 




e'''*s(a5(/))xmXi^0 j lis, cr)f{s, r/)/(s, - rj - a) dq da ds. 



The idea will be the following: on the support of the symbol of ()7.6bp (respectively (|7.6cp ). drjcj) 
(respectively da<j)) does not vanish, thus one should integrate by parts in r] (respectively a). On 
the support of the symbol appearing in ()7.6dp . both d^cp and dj^cp vanish, but the identity (j3.6p is 
the remedy: it essentially converts d^(j) into a combination of da(l> and drjcf), making the integration 
by parts in (r/, a) possible. 

Bound for ^7. 6b ) As explained above, we shall in order to treat this term integrate by parts in a 
usmg the identity t\t% ^'"^ = ^his gives (for the sake of simplicity in the notations, we do 
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not differentiate between standard product between scalars, and scalar product between vectors) 



i^XiXil^Q ( -77 ) d^ 'f{s, (y)f{s, ri)f{s, i - r] - a) dr] da ds 



{1.1.) W = IJje i^^^p — 



(7.7b) 
(7.7c) 




' di.(j)da4> 



^XlXifJ-Q ( ^ ) f{s, cr) f{s, 'n)da f{s, C - r] - a)dr]dads 



f{s, cr)f{s, r])f{s, - r] - a) dr] da ds. 



In order to estimate (j7.7ap . observe that tlie symbol 



'-XiXiP- belongs to 75/2- Indeed, the 



l^pyj singularity of d^cj) is canceled by /u, and the ^^_^^^^i/2 singularity of d^(j) is canceled by |g . 
Thus, applying first (IA.2j) - since the Fourier support of the integrand lies within a ball of radius 
comparable to - one gets 



Be 



ds 



8/5 



(7.8) 



< 



H\ie ||m|Ii6 ds 



< I s 4 s2'^o IIX/II2 ||ii||]^g ||n||-^g ds < ||u||j^ J s^s^^^^s^s »s s ds < 



where the last inequality holds provided that ^ + ^Sq + 6 — ^ < —1. Thus we choose 60 and 
such that this inequality holds, as well as (j7.4p . 

The estimate of (|7.7bp is almost identical; as for the estimate of (j7.7cp . it follows in a very similar 



way. Let us say a word about it: the symbol clo- 



the type A, 



-^;^XIXl^^ 



only contributes singularities of 



and 



treated; the same strategy is valid here. 



. We have seen many instances of how these singularities can be 



Bound for \7.6c\j It can be derived in an identical fashion. 



Bound for ( 7.6d) For this term, we use the identity (j3.6p . that we rewrite as follows: there exists 
a symbol p £ Tq such that 

(in the above, p is a matrix. It is linear in the arguments between brackets). 

Using the identity Xiiid^4' = j^Pi^iV-,^) [9r?e*'^'^, 5^6*^"^] to integrate by parts in (I7.6dp gives 
terms which can all be estimated as above. This concludes the proof of claim 17.11 

Bound for V^e~**^^''^(73 in L°° , < < 4. First, notice that we can not integrate by parts in 
time as we did in the weakly resonant case since now the phase vanishes on a large set. Actually, 
the proof of this bound follows the steps of the previous argument, in particular Claim [7T1 Indeed, 
in the previous argument, we derived a bound for ||x(73||2; by the same token, one might prove 
bounds on || V'^"'"^x(73||2. These two quantities barely miss to control 11^^(7311^3/2, which would 

give the desired L°° decay, by the dispersive estimate in Lemma lA.ll What is needed is a little 
bit of additional integrability, for instance replacing the Lebesgue index 2 by 2 — for a small 
constant 61. As we will see, this is possible if one is prepared to lose a small power of s; but a 
small enough power of s is harmless in that it does not prevent the integrals over s occurring above 
from converging. Hence, we will prove the same estimate as in the previous argument, but in L"^^^^ 
instead of L?. 
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In order to implement this program, we start by splitting into high and low frequencies: 



93 = J J e''<i'xif^e(^-^)f{s,rj)f{s,a)f{s,C-rj-a)dr^ 

r-t 



+ 



fis, v)f{s, cr)f{s, i-r]-(j)dr] 



def 



93,low ~l~ 93,high' 



Bound for V^e ^^"^^^^ 93 low , with < k < A. Rewrite gsiow as 



t 

93, low = I G3^iowis)ds, 

with the obvious definition for G^^iowis)- Due to the frequency localization of G^^iow, one might 
even write 



(7.9) 



93,low — j P<C So log sG 3, low {s) ds. 



Computing xG3^iow{s) gives terms similar to those appearing inside the time integral of (j7.2ap . 
dLlEl), (iT^Sbl) . namely 



(7.10a) HxG3,iow{s)) = 1 1 e'^^xwQ ( 4: ) f{s, <j)f{s, v)dif{s, ^ - v - a) drj da ds 



a 



So 



(7.10b) 
(7.10c) 



+ / / e''*d^ixifJ.)e ( ^ ) fis, a)fis, v)fis, i - - a) drj da ds 



a 



aSo 



+ 



j I e''ts{d^cP)xil^e (^) fis, a) fis, rj)fis,C-rj-a) dr^ da ds. 



It has been proved in the previous section that all of these terms can be bounded in L^. More 
precisely, an inspection of the proof there reveals that xG3^iow{s) can be written 



xGs^iowis) = e''''^'^^^ Gs^iowis) with G3^iow{s) 



< 



for some number ki > 0. By repeating the same argument (see for instance ()7.8p ). one can prove 
that, denoting by 61 a small enough constant whose precise value will be set later on, 



G3^lowis] 



< 



2-S, - .1+^ 



l7/l|3 

M\x- 



Lemma lA.ll gives 



V^P<C5o log sGs^ioiu (s) 



2-(5i 



V^P<C(5ologse*'^^ '' G3^iowis) 



(7.11) 



< f /o^2V-*i ^ 



G3^lowis) 



2-5i 
2-<5i 



\U\\xS " I s "s 



1 



S^ 2 
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Thus, choosing 5q and 5i small enough, and using successively the dispersive estimate in Lemma [A. H 
the equation (|7.9p . and the above bound. 



V e gsAow 



<i 

OO I 



1 



b\/^ ~ t 



(7.12) 



< 



||-f<C5ologsa;G3,/o«,||2_5j + 



2-(5i 



2-5i 



<7ll«llv 



1 1 

2-<5i 2 -L 



^ '^•5 £ -Il«ll V- 

1+^ ~ . II IIX 



Bound for V^e ^^^^'^ g^^high in L°° , with < A; < 4. Rewrite gs^high as 



gSyhigh 



i: 



e''^''"G^,Mgh{s)ds. 



Computing V^xG^^highis) gives terms which can be estimated just like ()7.3ap : one gets a bound of 
the type 



'^'^xGs^highis] 



< 



1 II ||3 



for a constant K2 > 0. Following the same arguments, but going to a smaller Lebesgue index 2 — 82 
gives, provided 62 is chosen small enough, and i < 6, 



"^^xGs^highis] 



1 



2-52 ^ s^+f '"^ 



3 

2 ll-llV- 



As above, this implies that 



V''G3,high{s] 



< _||^|||. 



p3/2 ~ 1 I ^ 
^1,1 2 



To deduce the L°° estimate, write 



V e gsMgh 



V^e'i^-'^^'''G^,h^gh{s) ds+ [ V^e'^^-'^^'^'G^^Mghis) ds. 

Jt-1 



The second summand of the above right-hand side can be estimated using Sobolev embedding and 
proceeding as in (j6.5p . As for the first summand. 



(7.13) 



y^e^^'-''^^''"G^,Mgh{s)ds 



< Mx 



t-i 



1 1 



ds < -||u||^. 



8. Estimates for terms of order 4 and higher 

In this section, we derive the desired estimates in L°° and Lp'{x^dx) for terms of the form (j4.3ap 
and (|4.3bp . These terms gather the contributions of order 4 and higher (in u) to the nonlinearity, 
once the normal form transform of Section U] has been performed. 

Being of high order, they decay very fast and can be estimated by brute force, leaving aside the 
question of resonances. Since a straightforward approach is sufficient, we only illustrate it in the 
case of ()4.3ap . which we denote 

g^= / e''^'^'R{s)ds. 
J2 
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Bound for V'^e'^'''^^'^ in L°°, for < A; < 4. We split 



R{s)ds 



into J2* ^ + //_!• The first summand is bounded via the dispersive estimate of Lemma lA.ll and 
Proposition IF.ll 



t-i 



V'^e'^'-'^^ ' R{s)ds 



< f ' -^\\R{s)\\ .ds<-\\u\\\. 



The second summand is bounded using Sobolev embedding and Proposition IF.ll 

rt 



V^e 



R{s)ds 



t-i 



< 



t-1 



R{s) ds< / \\R{s)\\h2+, ds<-\\u\\\. 

t-i ^ 



Bound for xg4 in L . The L norm of xg^ can be bounded as follows: 

ft 



\xg4.\\2 



< 



X I e""^^^^ Ris) ds 



\^^e-^'''R{s)ds 
t rt 



+ 



e^"^'^'xi?(s) ds 



^ J2 '^"^^'^'*"4/3 + WxRW^ds. 



The first term in the last line above can be estimated directly using point [i) of Proposition IF.ll 
we leave this to the reader. The bound for the second term follows with the help of point {ii) of 
Proposition IF.ll 



t ft rt 

\\xR\\2ds< / \\{x)u\\2\\u\\l^,■i,oo ds < 

2 J2 J2 



\u\\2 + S 



+ \\xf\\, 



I'ull^a.oo ds < 



9. Scattering 

In this section, we prove the scattering result (Corollarv ll.ip . On the one hand, decomposing / 
into terms of the form (j4.H44.3bj) . and examining them separately, it appears that it is a Cauchy 
sequence in in the sense that there exists a constant k > such that 

\is<t<2s, ||/(t)-/(5)||2<l. 

On the other hand, we know that 

ll/WII 

Interpolating between these two inequalities gives, for a constant (7, 

iis<t<2s, \\f{t) - f{s)\\HN-Nc,s^LH.^-c,s,,) < 
Choosing Cq properly makes f{t) Cauchy in time, which gives the desired result. 
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Appendix A. A few results of linear harmonic analysis 
Littlewood-Paley theory. Consider 9 a function supported in the annulus C(0 
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, |, |) such that 



Also define 



Define then the Fourier multiphers 

2^ 



^3 ~ 



p>i = 1 - e 



and similarly P<j, P>j- This gives a homogeneous and an inhomogeneous decomposition of the 
identity (for instance, in L?) 

^Pj=ld and P<o + ^Pj = ld. 

All these operators are bounded on spaces: 

ifl<P<oo, ||P,-/||p<||/||p , ||P<,-/||p < ll/llp and ||P>,-/||p<||/||p. 
Furthermore, for Pjf, taking a derivative is essentially equivalent to multiplying by 2^: 

if 1 < p < GO and a e M, P"-P,7||p ~ 2"^' ||Pj/||p 
if 1 < p < GO and £ e Z, \\V^Pjf\\p ~ 2^^\\Pjf\\p. 
Also, we recall Bernstein's lemma: ifl<q'<p<oo, 

(A.2) \\Pjfl < 2'^'(^i) and ||P<,/||^ < 2'^(^^) ||P<,-/||^ . 

Finally, we will need the Littlewood-Paley square and maximal function estimates 

Theorem A.l. (i) If f = J2fj> '^^^h Supp(/j) C C(0, c2~-^ , C2~-') (the latter denoting the annulus 
of center 0, inner radius c2~^ , outer radius C2~^ ), and \ <p < oo, 



(A.l) 









1/2 




< 






3 


p 


3 





Furthermore, denoting Sf 



def 



Y.^P3f? 



1/2 



, ll^/llp~ll/llp- 



(ii) Ifl<p<oo, denoting Mf{x) =^ sup \P<jf{x)\, \\Mf\\p < 



Fractional integration and dispersion. 

Lemma A.l. (1) (Fractional integration) If 1 < p,q < oo, < a < |, and ^ — p ^ then 
1 



f 



< 



(2) (Dispersive estimate) 



< 



1 



oo 
2 



i}3/2- 

1,1 

2 2 



(3) // 1 < p < 2 < g < oo, < a < ^, and a = ^- ^, then 



< 



p- 
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(4) Ifl<p<2, £>0 and 2H^ > I, 
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f 



< 



p 



2^0 {2H'') U 



1 i_ ' 

2 , 



Proof. We only prove the second and forth points, the rest being standard or elementary. 
Proof of (2). The stationary phase lemma gives 



<1 
~ t' 



This implies that 



By scaling, 



< 



l^o/lli 



1 



m <tII^'o/iIi 

OO b 



OO I 



This inequality gives immediately the desired conclusion. 

Proof of (4). First notice that it suffices to show this result for ^ = 0. It will follow from interpolation 
between the estimate, which is clear, and the estimate, which reads 



if > 1, 



By scaling, it suffices to prove this estimate if t = 1 and j > 0. This is done as follows 



< 



1/2 



k=l 

j 



1 + E 



k=l 
j 



2k 



k=l 



1/2 



1/2 
2 

1/2 



1/2 
2 



2k 



1 ] ^W'/' 
2k ' ^ 



1/2 
2 



1/2 



1/2 
2 



fc=i 



□ 



Appendix B. Some general facts on pseudo-product operators 

Let us first give the definition of pseudo-product operators, which were introduced by Coifman 
and Meyer; we only consider the bilinear and trilinear cases, which are of interest for our prob- 
lem. Bilinear (respectively trilinear) operators are defined via their symbol 771(^,77) (respectively 
'm{i,ri,a)) by 

Bm(^,n,a)ifi, f2, fs) J^'^ j m{i,r],a)fi{a)f2{'n)h{i-ri-(^)dvda. 

The fundamental theorem of Coifman and Meyer |CM78] states, under a natural condition, that 
these operators have the same boundedness properties as the ones given by Holder's inequality for 
the standard product. 
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Theorem B.l. Let n be either 2 or 3, and suppose that m satisfies 

1 



(B.l) . . . 9?"m(ei, ...,in)< 

for sufficiently many multi-indices. Then 

\\Bm{fi,-- • ,/n)|lr < ll/illpi • • • ||/n||p„ 7- = ^ ^ 1 < Pi . . . Pn < oo and < r < CO. 

The class of symbols allowed by the Coifman-Meyer theorem (typically, symbols homogeneous 
of degree and smooth outside the origin) does not contain the symbols occurring in our analysis 
of the water wave problem. We will therefore be led to introducing a new class of symbols, in 
section O 

In the trilinear case, these new symbols will amongst other discrepancies with the Coifman-Meyer 
case exhibit flag singularities, to which we now turn. What is a flag singularity? If one considers 
two bilinear symbols mi and m2, and one trilinear symbol m^, all of Coifman-Meyer type, the new 
symbol m{E,,r],a) = ?Ti3(^, r;, cr)mi (^, 77)7712 (^ — rj,a) does not satisfy the Coifman-Meyer estimates 
(in a neighborhood of ^ = 7/ = and ^ — r/ — (T = cr = 0). This type of symbol is said to have 
a flag singularity. The following theorem gives the boundedness of such operators. A result of 
Muscalu |MU07j is slightly less general, but we are able to give a simpler proof since the range of 
Lebesgue exponents we are interested in is smaller. 

Theorem B.2. 7/777-1(^,77), m2{i,rj), m3{£,,rj,a) are symbols satisfying the estimate iB.l]) . then 

(/,5,/i)|| <II/IUII 

for 1 < p,pi,P2,P3 < °o ^^'^ ^ + ^ + ^ = |- 

Proof. Step 1: Partition of the frequency space. Set m{S,,r],a) = ms{£,,r],a)mi{^,ri)m2{(, — ri,a). 
First observe that there are certain regions of the (^, 77, a) plane where 777 satisfies the Coifman- 
Meyer estimates (jB.ip : then the Coifman-Meyer theorem applies, and the desired estimate is proved. 
Thus, using a cut-off function, we can reduce the problem to the regions where the Coifman-Meyer 
estimate (jB.ip does not hold for 777, namely 

Ai U A2 =^ m + |r?| « k|} U {|C - 77I + \a\ « m, 

We further observe that on Ai (respectively: on A2), m2{C — rj,a) (respectively mi{^,r])) satisfies 
the Coifman-Meyer estimate in {^,r],a). Thus with the help of cutoff functions, we can reduce 
matters to symbols of one of the two following types 

(B.2a) XAi{^,r],cr)m3{£^,r],a)mi{^,r]) (where XAi localizes near Ai) 

(B.2b) XA2(C) o-)"T'3 (■?,??, o")7n.2(^ - ??, o") (where XA2 localizes near A2) 

where 7713 stands for m^m2 or 77737771. Estimates for one of these symbols can be deduced from 
the other by duality and using that for any two symbols /7 and and (•, •) denoting the standard 
(complex) scalar product, 

{Bli(i,r),a)u{i,ri){fl^ f2, fi) , fA) = {B^(^^tj,(,-'q-a-E,)'^(~cj,i-r]-a){fA,fi,f2), fl)- 

Therefore, we simply prove estimates for the symbol ()B.2ap . We define the cut-off function XA-i by 
(B.3) B.^^^(^^^^„){f,g,h)= ^ P^k^iooiPkfPk'h) P<k^ioog. 

\k-k'\<l 

We will suppress the index k' in the sequel and just take k' = k to make notations lighter. 
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Step 2: series expansion for the symbol Expanding the symbol 171^(^,7], a) near {rj,(,) = gives 

M-l 

\a+l3\=0 

where ^a,i3{(^) = ^|g^"i3(0, 0, a) has homogeneous bounds of degree —\a + /3| i.e. satisfies 
|5j<l*Q,/3| < |(7|~l'>'l~l"+^l, and the remainder R is such that 

(B.4, |<,a,.«,„.,|.o(M±l^). 

This estimate implies that, if M is chosen big enough, the symbol resulting from the multiplication 
of R and miXAi is of Coifman-Meyer type, thus we can forget about it. Thus it suffices to treat 
the summands of the first term of the above right-hand side; to simplify notations a little in the 
following, we simply consider the case a = 0, and therefore get symbols of the type 

(B.5) m(^,r/,cj) = XAi(^,r/,o-)^>i(cr)r?*mi(^,r?) , 

where has homogeneous bounds of degree —i, and mi satisfies the Coifman-Meyer estimates. 



Step 3: Paraproduct decomposition for the symbol mi. Proceeding as in the original proof of 
Coifman and Meyer |CM78j . we will now perform a paraproduct decomposition of m, and then 
expand the resulting symbols scale by scale. Write 

Brn, (/, g) = Yl ^-1 iPjfP<J~l9) + Yl ^-1 {P<J-lfPj9) + (^J-^' ^^9)- 

3 i li-^l<i 

Next consider the symbol of one of the elementary bilinear operators above, for instance -B^i (Pj': P<j-i')- 
Denote m\(^,r]) for this symbol, which is compactly supported (in (^, and expand it in Fourier 
series 

where we denoted c for a constant, x fo'^ ^ cut-off function, and Op^q for the Fourier coefficients. 

It is now possible to forget about the summation over p, q because of the fast decay of the 
coefficients Op^q, which results from the smoothness of the symbol. Indeed, the complex exponentials 
g«c2 become in physical space translations, which add polynomial factors to the estimates 

to come. These polynomial factors are offset by the decay of the ap^g. This leads to replacing ap^q 
by aj S It is now possible to consider that the aj are actually constant in j: if this is not the 
case, it essentially corresponds to composing one of the argument functions by a Calderon-Zygmund 
operator bounded on Lebesgue spaces, which is harmless. For these reasons, it will suffice to treat 
the case when mi corresponds to one of the three paraproduct operators 

if, 9)^ YP,fP<,-i9 ; j;p<,-i/P,5 ; Y.PifPi9 
j j j 

(we suppressed the index £ in the last summation to make notations lighter by taking i = j). 
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Step 4- Derivation of the model operators. Combining this last line with (|B.5p . we see that the 
operators of interest for us become 

J2P<k~woPj {MD)PkfPkh) V^P<,_iP<fc_iooff 
j,k 

g-) ^ -P<fc-ioo-P<i-i i^i{D)PkfPkh) V*PjP<fe_iooff 

Y,P<k-iooP, iMP')PkfPkh) V'PjP^k^ioog. 

We now make some observations which allow us to simplify the above operators: 

• First remark that ^i{D)Pj^, ^^Pj ^^'^ ^^P<j '^^^ written respectively 2~'^^Pk, 2^^ Pj and 
2*-'P<j with obvious notations. Since the operators with tildes have similar properties to 
the operators without tildes, we will in the following forget about the tildes. 

• Next notice that due to the Fourier space support properties of the different terms above, 
it is possible to restrict the summation to j < k — 97. 

• Finally, since P^k-iooPj = Pj and -P<j-i-P<fe-ioo = P<j-i for j < k — 103, it is harmless to 
forget about the P^k-ioo operators in the above sums. 

All these remarks lead to the following simplified versions of the above operators: 

^ 2^(j->^)p^(P,fP,h)P^,^,g 

k>j +97 

(B.7) E '^'^''^^P<o-i{PkfPkh)Pjg 

k>j +97 

2'^^-^^P,{PufPkh)P,g 

k>j +97 



Step 5: The case i = 0. If i = 0, observe that, due to the Fourier support properties of the 
Littlewood-Paley operators, the operators in ()B.7p are equal respectively to 

(B.8) ^P<,_i(^^P,/P,/i^P,-5 



up to a difference term which is Coifman- Meyer. But the operators in (jB.SP are simply compositions 
of bilinear Coifman-Meyer operators. Thus the desired bounds follow for them. 
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Step 6: The case i > 0. If i > 0, we see that it suffices to prove uniform estimates in J > for the 
operators 



(B.9a) 
(B.9b) 
(B.9c) 



3 

Y,P<3-iiPj+jfPj+jh)Pjg 

3 

Y,PjiP3+jfP3+jh)P,9 



since the desired result fohows then upon summation over J. The estimate rehes on the Littlewood- 
Paley square and maximal function estimates (Theorem lA.ip . and on the vector valued maximal 



function estimate 
for (iRQi]) 



1/2 



1/2 



~ II (Z^i/i^j lip (s6e Stein [S93j . chapter II). This gives 



dEliDllp 



Y,[P3iP3+jfP3+jh)P<,-l9Y 



< 



\\Mg\\ 



P2 



1/2 



1/2 



< 



1/2 



Mg [P3 {P,+jfP,+M 



Y,[P3iP3+jfP3+jh)Y 



^ Ibl 



P2 



PP2 
P-P2 



1/2 



Y,[M{P,+jfP,+jh)f 



^ II5I 



P2 



1/2 



E [P3+jfP3+M 



^ II5I 



P2 



1/2 



Mf [P3+jhY 



^ \\9\\p2 IIM/IU 

The terms (iBlgbjl and (iRQcj) can be estimated similarly. 

Appendix C. Analysis of a class of bilinear pseudo-product operators 



□ 



In this section, we define new classes of bilinear pseudo-product operators, which occur in the 
analysis of the water wave problem; then, we prove boundedness of these operators. 

Definition of the classes Bg and Bs- Before defining these classes, it will be convenient to think 
of a symbol in a more symmetric way than we have been doing so far. The motivation is the 
following: in order to prove a bound for the bilinear operator Bm{f,g) in U one uses a duality 
argument and proves that for h £ L'^ , we have / Bm{f, g)hdx < 00; this shows that the Fourier 
variables ?7, ^ — — play symmetric roles. 

def 

Given a symbol m{^,rj), we can as well write it rn^^jS, — 1]) — '"^(C)^) 01 more generally as 
a function of two of the three variables (r/, ^ — 'r],—^^). By an abuse of notation, we will denote 
indifferently m for all these symbols. Thus, denoting {(,1,^,2, £,3) for the three Fourier variables 
(rj, £ — f], — ^), and picking two indices i and j, we can always write m = m{^i,^j). 

Definition C.l. A symbol m belongs to the class Bg if 

• It is homogeneous of degree s. 

• It is smooth outside of {^1 = 0} U {^2 = 0} U {^3 = 0}. 

• For i = 1,2,3, if << |Ci+i|, |^i+2| ~ 1; it is possible to write m = A MCil^^^, ]|^,^i+i) • 
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(where we use the convention that ^4 = ^1 and ^5 = ^2^- 

Thus, roughly speaking, symbols in Bq are of Coifman-Meyer type except along the coordinate 
axes, where they are allowed to have a singularity like Mihlin-Hormander (linear!) multipliers. 

Notice that given the boundedness properties of Mihlin-Hormander multipliers and Coifman- 
Meyer operators, Theorem lC.il of next section, which gives boundedness of bilinear operators with 
symbols in Bq, should be no surprise. 

We now define a new class of symbols, which corresponds to bilinear operators of paraproduct 
type. 

Definition C.2. A symbol m belongs to the class Bg if 

• It belongs to Bg- 

• It satisfies the following support property: Suppm(^,ry) C > |^|}. 

The interest of the class Bg is the following: taking derivatives of B^{f,g) corresponds to mul- 
tiplying it by ^ in Fourier space, li m G Bg then the support restriction on m means that ^ is 
always dominated by r]. Thus one expects (see next section for a precise formulation) something 
like |V'=i?„(/,<7)| < iV^fM- 

Calculus with symbols in Bg and Bg. We begin with the action of derivatives on Bg. 

Lemma C.3. (i) If £ Bg, one can write 

1 1 

d^fJ'iv, = l^^ + TTil^^ + 77 tA*^ with (/i^ , /i^) G Bg-i X BgX Bg 

(a) If fj, £ Bg and v £ Bg' , then fiu G Bgg' . 

Proof, (ii) follows from the definition of Bg. To prove (i) one notes that d^fi contributes /xi if |r/| 
is the smallest variable, it contributes |^ if I'?! is the smallest variable, and it contributes j^z^ if 
— r/l is the smallest variable. □ 

Finally, the following theorem gives the crucial boundedness properties of operators with symbols 
in Bg and Bg. 

Theorem C.l. (i) If m belongs to the class Bq, 

\\BUf,9)\\r ^ ll/llpll^llg ^/^ + I = ? and I <p,q,r < 00 

(ii) If m belongs to the class Bg, and if k is an integer, 

V''Bm{f,g)\l < \\A'^'fUg\\g if 1 + 1 = 1 andl<p,q,r<^ 

Remark. Since Riesz transforms are not bounded on L°°, the statement of the above theorem 
becomes wrong if any of the indices p,q,r is 00. We need to go around this difficulty when deriving 
estimates for the water wave system, and this unfortunately makes estimates a bit longer. 

Proof. To begin with (i), let m be a symbol in Bq. Away from {^1 = 0} U {^2 = 0} U {^3 = 0}, 
the Coifman-Meyer theorem [CM78j can be applied to m and gives the desired result. We will 
simply consider the case |r/| << |^| (the other cases can be reduced to this one by duality), thus 
we consider in the following 



B'if, g) ^= Yl Bm (P<i-ioo/, Pjg) . 
3 
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By definition of the class Bq, m can be written in the set \'r]\ « |^| as m(S^,rj) = A 

or, by homogeneity, r]) = A (j^\T/2^ • Expanding this expression in j^j^^^ yields 



11^ \fj\k/2 / ^ \ 



First notice that the m/; are smooth and homogeneous of degree 0. Second, observe that for M 
big enough, the singularity of the remainder at ry = becomes so weak that the remainder satisfies 
estimates of Coifman-Meyer type; thus we take M big enough and forget about the remainder. 
Expanding in spherical harmonics (which we denote (Z^)^gj^) leads to 



AI 



m 



By the Mihlin-Hormander multiplier theorem, the operators Z£ (^) are bounded on Lebesgue 
spaces, with bounds growing polynomially in i; on the other hand, the smoothness of m entails 
fast decay in of the coefficients ak//'- In the end, we can thus ignore the summation over £, 

i' and k, and matters reduce to 



It is then routine to get the desired estimate using slight extensions of the Littlewood-Paley square 
and maximal function estimates (Theorem lA.ip 



(C.l) 



-k/2 



P, 



<j-100 



< 



1 

2\ 2 



< 



< 



sup,- 



2-"-l'-P.,i.,„A>-l^f(x)\(j^iP,gf 



1/2 



sup,2-^'=/2p<,_ioo^'/V 



1/2 



< 



In order to prove (ii), consider a symbol m in iS^, and simply observe that 

Due to the support condition satisfied by m, the symbol |^p.^^ m belongs to B , therefore applying 
(ii) gives 



V'=S,n(/,5) 



< 



□ 



Appendix D. Analysis of a class of trilinear pseudo-product operators 



In this section, we turn to the trilinear operators which occur in the analysis of the water-wave 
problem: we define new classes of symbols adapted to them, and prove their boundedness. 
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Definition of the classes 7^ and Tg. As in the quadratic case, it will be convenient to put on 
an equal footing the four Fourier variables S,,rj,cr,(, — r] — a. We adopt the following convention: 

(1) First, we call oi = — ^, a2 = r], = a, = ^ — t] — a. 

(2) Then, we partition the rj, a) space into regions where the (|aj|) are essentially ordered (in 
other words, for each of these regions there is a permutation a such that |ao-(i)| < |ao-(2)| < 

|a<T(3)l ^ |acr(4)l- 

(3) Finally, we set = a^fj^ . 

In other words, (^i) is a convenient labeling of — ^, r], a, ^ — rj — a since it satisfies |^i| < |^2| ^ 
161 < 1^41 and 6 + 6 + 6 + ^4 = . 

As in the quadratic case, we abuse notations by denoting indifferently m for the symbol m(^, rj, a) 
or its expression in any coordinate system, for instance m'(6)6i6)- 

Definition D.l. A symbol m belongs to the class Tg if 

• m is homogeneous of degree s. 

• m is smooth outside a conical neighborhood Oij D {^i = 0} U {^j + = 0} for {i,j) = 
(1,2), (1,3) or (2,3). 

• For\Ci\ ~ 161 « 161 ~ 1, and |6 + 61 ~ |6| if (iJ) = (1,2), m = ^ (6, 6) ^ (6, 6, 6) 
(flag singularity). 

. For 161 « 161,161,161 "1 = ^(161'/', 1^,6,6). 

• For\Ci+Cj\ « 1^11,161,161,161 m = ^(|ei + 6f/',|T|^, 6,6,6) ■ 

. For 161 « 161 « 161, 161 -i,m = A 1^, |6l^/^ 1^, 6) • 

• For 16 + 61 « 161 « 161, 161 ~ 1, m = ^ lirfc' l^il'^'' l^'^s) . 

Remark. One should think of symbols in Tq as being of Coifman-Meyer type, except that they 
might exhibit flag singularities, and they are allowed to have singularities like Mihlin-Hormander 
(linear!) multipliers along the coordinate axes, and along one axis corresponding to the sum of two 
coordinates being zero. 

Though this definition is quite involved, it somehow corresponds to the simplest class containing 
all the symbols which occur in the analysis of the water waves problem. For instance, symbols of 
the type ~ V?^) (which appear in Section 6) contribute singularities of the type ||| if 

one of the coordinates vanishes, flag singularities, and singularities along one axis 6 + 6 • This 
already accounts for nearly all points of the above definition. 

How does one prove that these symbols are associated to bounded operators (Theorem \D.1\) ? It 
is important in the proof that symbols in Tq can, by power or Fourier series expansions, be reduced 
to tensorial products of functions of one 6 only; in other words it is possible to separate variables. 
Notice that symbols of the form (for instance) could not be treated by such a method; but such 
a behavior is not possible in the class Tg. 

We next define a new class of symbols, which somehow corresponds to paraproduct operators. 

Definition D.2. A symbol m belongs to the class Tg if 

• It belongs to Tg. 

• It satisfies the following support property: Suppm(6 ??, o") C {\a\ > |^|, \r]\}. 
Calculus with symbols in Tg and Tg. We begin with the action of derivatives on Tg. 
Lemma D.3. (i) If fi £ Tg, one can write 

a^/x(6 V, = /^^ + TTT/^^ + T7-—il^'^ + \ l^'^ + 77 " 1/^^ 

1^1 l^-r/l \i-CF\ 
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with (^i,;u2,/x3,/i^,;u^) G Ts^i xTsXTsXTsX 

(ii) If G Tg and v ^Tgi, then ^lu E Tggi. 

Proof. Follows from the definition of Tg. Actually, due to the fact that /i is smooth outside a 
neighborhood of Oij for (i, j) = (1, 2), (1, 3) or (2, 3), we also have that /Us = or /i4 = 0. □ 

Finally, the following theorem gives the crucial boundedness properties of operators with symbols 
in Tg and Tg. 

Theorem D.l. (i) If m belongs to the class Tq, 

\\Bniifl, f2, f3)\\r ^ II /l llpi II /2 ||p2 II /s Ilp3 ^/ FT + ^ + ^ = r "-"^^ 1 <Pl,P2,P3,r < OO 
||^m(/l,/2,/3)|lB0_^ ^ Il/l||pill/2||p2ll/3||p3 ^/ FT + ^ + ^ ^ 

(ii) If ni belongs to the class Tg, and if k is an integer, 

V'^i?™(/l,/2,/3) < P'+7l||pill/2|UI|/3lU ifv,+V,+V,=-r (^ndl<pi,p2,P3,r<00 



r 



V'^i?™(/i,/2,/3) < P'+7ilUII/2||p2ll/3lU ^/^ + i; + ^ = l andl<pi,p2,P3<oo. 



Remark. The estimates above involving Besov spaces come as substitute for estimates in , which 
are wrong (for instance, the Riesz transform is not bounded on L^). One can get in the same way 
substitutes for estimates in L°° . 

Proof. The proof of (i) is very similar to the one of Theorem IB.2I and Theorem lC.il first, partition 
the frequency space in order to distinguish the different regions appearing in the definition of Tg. 
Then, in each of these regions, expand the symbol in order to separate variables; finally, perform 
the desired estimates using the Littlewood-Paley theorem. 

In the region ~ |^2| << |?3| ~ 1, and |^i + ^2| ~ |Ci| if ihj) = (li^), it suffices to apply 
Theorem IB . 2 1 on symbols with flag singularities. 

Let us also sketch how one deals with the region << |^2| << iCsl, l^il- Expanding the symbol 

('^^pT' iIti) l'^2r^^) il^yjCs^ ill power series in j^^j^/^ and |'^2|^^^; and in spherical harmonics in 

k/2 



il^ and il^ gives 

l?l| l?2| " 



where the functions ^kk'ee' are homogeneous of degree — -y and decay fast with the indices by 
smoothness of A. 

The above sum over k, k' can be taken to be finite, the remainder giving a Coifman-Meyer 

/ 1^ I \ k/2 j^, 12 

operator; furthermore, it is easily checked that the powers ( ^ j and |^|2 appearing above 

cancel with the homogeneity of ^kk'ei'] thus it is possible to ignore the sum over k and k' . 

As for the sum over we rely on the fast decay of the functions ^kk'U'i which offsets the 
polynomial growth of the bounds (in Lebesgue spaces) of the Fourier multipliers Zi and 

Zc ^]|^^ ; thus it is possible to ignore the sum over i and £' . 

Finally, up to rotation of the Fourier variables, and a duality argument, it is possible to assume 
that (6,6,6) = (<7, ^ - - cr). 

The above considerations lead to the model operator 

B{f,g,h)= Yl P<f -100 fPf 9 Pjh. 

0<k,k'<N j j'<j-100 
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The boundedness of this operator is easily estabhshed: it is essentially a composition of paraprod- 
ucts. 

The extension to Besov spaces based on follows from boundedness from to Bq ^ of Mihlin- 
Hormander type Fourier multipliers. 

Finally, the point (n) is proved just like for Theorem lC.il □ 

Finally, we need the following proposition, which combines fractional integration and flag singu- 
larity: 

Proposition D.4. (i) If m{S,, rj^a) £%, < a < 2, and 1 < pi,p2,P3,r < oo, then 

5^_L_^„(5,^,^)(/l,/2,/3) ^ Il/l||pill/2||p2ll/3||p3 /^^^^ FT + ^ + ^ ~ ^ = 7- 
(ii) If m{S,,rj,a) G %, < a < 2, and 1 < pi,p2,P3,r < oo, then 



V'=i?_^„(^,,,,)(/l,/2,/3) ^<P'+7llUII/2||p.||/3||p3 for^ + j-,+ 



2 L = 1 

P3 2a r ' 



Proof. The proof follows the pattern of the proofs of Theorem IB .21 Theorem lC.il and Theorem lD.il 
In the end, it thus somehow reduces to the case m = 1, for which the estimate is clear since 

B\^-r]\~^'{fl: f2, fs) = /2^(/l/3) 



and thus, by Holder's inequality and Lemma lA.l| 

1 



B 



1 

li-vi' 



.(/i,/2,/3) <||/: 



Ip2 



(/1/3) 



< 



P2- 
P2r 



II/iIIpJIMIpJI/sI 



P3- 



□ 



Appendix E. Proof of proposition 12.21 

Let a = {x, z) £ and b = {y,h{y)) G S. Then ipT-i can be represented by the double layer 
potential 

V'w(a) = J Kb)N ■ VG{a - b)dSib) = ^/x(6o) + J (Kb) - Kbo))N ■ VG{a - 

s s 

where G{a — b) = — is the Newtonian potential and 60 is an arbitrary point on S. By 

defining K{x,y) = ^/ 1 + \V h{y)\'^ N {b) ■ VG{bo - b) for bo = {x,h{x)) € S and 6 = G S 

then 



(E.l) 



' Kix y) = 

i-„(\x-y\^ + {h{x)-h{y))-^)i 

/ \K{x,y)\{l + \x- y\^)dg < \\h\\w2,^ + ||V/i||lp < Eq 



for 2 < p < 4 and from standard singular integral calculations |F095| as a — > 60 
(E.2) ^ 



^fJ.{x) + J fi{y)K{x,y)dy = ipix). 
From (|E.ip it is clear that K maps 

^ Loc j^j^^ ^ for < a < i. Therefore we can 
solve for /i by a Neumann series expansion to obtain 



(E.3) 
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To obtain estimates on derivatives of /i, we change variables from y to z = x — y in (|E.2p and write 
J{x,z) = K{x,x + z). Since dxJ{x,z) satisfies inequality (IE.1|) with one additional derivative on 
h, i.e., 

J \d^J{x,z)\{l + \z\^)dz < ||/i||h^3,oo + WVhWw^.r, < Eo, 

then differentiating ()E.2p with respect to x gives 

-dfM{x) + f dfj,{x — z)J{x, z)dz + f ^{x — z)dxJ{x, z)dz = d'ip{x) 



which implies ||9;u||j;,oo < H^-^Hioo + ||^|| . i < HS^/'Ulcx) + ||'0|| . i by (|E.3p . Repeating the above 
argument twice we obtain 

wwc^ < \mc2 + 



■ 1 . 

C3 



Next we estimate Mtl^. To do this we fix a point bQ £ S and use normal coordinates in a neighbor- 
hood of S to restrict a near the boundary to the line a = + vN{hQ). Thus 

N{bo) ■ VVw(a) = J ifiib) - fiibo))D^GiN{b),N{bo))ia - b)dS{b). 
s 

For |6 — 6o| large and u small \D'^G{a — b)\ < |6 — 6o|~^ and thus the above integral can be bounded 

by \\fi\\ ■ 1 . For \b — bo\ small we write 

c ^ 

N{bo) = eib,bo)N{b)+j{b,bo)T where r G T^^S. 

The term involving r is integrable due to the vanishing of 7(60! ^o)- By repeating the argument 
that led to inequality (|E.3P we obtain 

1 1 {i,{b) - i,{bo))D^G{N{b),jT){a -b)dS{b)\ < My/2 < IIVIIci/- 
s 

The kernel of the remaining term 



1 = J {fi{b) - ^i{bo))D'G{N{b), N{b)){vN{bo) + 60 - h)dS{b) 
s 

is hypersingular as u ^ Q and can be dealt with by using the identity 

= AG = AcjG + K,N ■VG + D'^G{N, N) 
for v < 0. This allows us to re-express / as 

1 = J Vn{b)VG{uN{bo) + 60 - 6) - in{b) - n{bo))KN{b) ■ VG{uN{bo) + bo-b) dS 
s 

which is a singular kernel and can be bounded as before 

W<M^i+IHI^i <||avil^i+M-. 

By repeating the above argument after applying tangential derivatives to N{bQ) ■ Vip we obtain 
This proves inequality (j2.ip . 
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Appendix F. Estimates on the remainder term R 

Recall that R is the remainder term defined at the beginning of Section [31 More explicitly, 
recalling that u = h + iA^^'^ip, R consists of the terms of order 4 and higher in u of the nonlinearity 



(F.l) 



G{h)ip + iA^/"^ 



2' ■ ' 2(1 + \Vh\ 
The following proposition gives bounds on R: 

Proposition F.l. // ||n||(4/3,oo is small enough, 



(F.2) 



i) 


V'^R 


< 












r 









u 



14/3,00 



I + I = r and I <p,q,r < 00 



(ii) \\xR\ 



2< 



\X)U\\2\\U\\y^i,. 



Remark. Both of the above estimates say that R behaves roughly like a four-fold product of some 
derivatives of u. Both of them are far from being optimal, except for one point: in estimate {ii), 
the weighted norm on the right-hand side does not carry any derivatives. 

We will only prove {ii), the estimate (i) being much easier. Furthermore, we observe that if the 
bound {ii) can be proved with R replaced by the fourth-order terms of G{h)ip, it follows easily for 
R. Indeed, a look at the nonlinearity (jF.ip shows that R essentially consists of products of G{h) 
with derivatives of u, but this is easily treated. Thus we will in the following prove {ii) with R 
replaced by terms of order four and higher of G{h). 



(F.3) 
Then 
(F.4) 



Given a harmonic function ip-j-i^ in J7, it can be represented via a single layer potential p a^ 

1 



il^nix.z) 



27r 



p{y){\x - + \z- h{y)\'')-^l^ dy if (x, z) G n. 



G{h)ij = p 



1 

2^ 



Vh{x) ■ {x-y) + h{y) - h{x) 
Pyy)- 10 . I, / N , . dy 



(|x-2/|2 + |/i(x)-/i(y)|2)3/2 

Notice that the kernel is very similar to the one given in (|E.ip 

Expanding (|F.3P in h, we get after evaluating on the boundary {x,h{x)) and applying A 

\h{x) - h{y)\^^ 



A^{x) = p{x) + ^ UnA J p{y) 

n>l 

'iip + Y^K^{p) 



\x-y 



2n+l 



dy 



(recall that yl = yl). Using Neumann series, p can be expressed as 

p = A'4)+ ^ f5mi...mNKmi ■ ■ -Kjni^Atp 

AfGN,(mi...miv)GN^ 



where \(3mi...mN\ ^ C"^^^ hmjv_ jj^ order to obtain weighted estimates, we will need the following 
claim, whose proof we postpone for the moment. 



Claim F.2. Ifn>l, 

(F.5a) 

(F.5b) 



\xKnAi;\\2<C''\\{x)u\\^\\uf^^X 
\{x)K^Zh<Cmx)Zh\\hf^^X 



^Note that in our representation p differs from the standard p by a factor of + |V/ip. 
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Assuming the claim, Proposition IF. II follows. Indeed, the contribution of in (|F.4p to R is 
given by 

^ ^ Pmi ...rriff -^mi • • • -^m ffAljj, 

miH \-mpf>2 

and this can be estimated using the claim. In order to estimate the contribution of the second 
summand of ()F.4p . the same procedure can be applied: expand the kernel in powers of h, and use 
the bound on the expansion of p in powers of h which follows from the claim. 

Proof of flF.5ah in Claim lF.21 Step 1: A paraproduct inequality. 
Lemma F.3. We have for A> and B > 1, 

(F.6) ||xV^yl/V^<7ll2 < \\A'/^f\\wA+B+i.^ \\{x)gh + \\{x)Ay^fh 

(F.7) ||xV^+VV^<7ll2 < ||/||vKA+s+2.oo \\{x)gh + \\{x)fh \\9\\w^+b+2,o. 

Proof. We focus on the first inequality, the second one being proved in an identical way. Using the 
paraproduct decomposition, we write 

xV^AfV^'g = x[Y^ V^APjfV^P<,g + V^AP^^fV^P^g 
j j 

Notice that 

xV^AP^.f = A^'"2P^jA'/^f + V^A'^^P^^xA'/^f 

(where we denote P<j for a symbol of the type A (^), with A smooth and supported in -6(0, C)). 
Thus 

xV^AfV^'g = ^ V^APjf{V''-^P<jg+V''P<jxg)+{A^-'2P^jA^/^f+V^A^/^P<jxA^/^f)V''Pjg. 
j 

This implies 

||xV^yl/V^5ll2 < 5;2^-(^+i)||P,/|U(2^-(^-^)||5||2 + 2^-^||x<7||2) 
j 

+ (2^'(^-i/2) |Ui/2/|l +2i(^+5)||xyli/2/||2)2^'^||P,5||oo 
< \\A'/'f\\w^+B+^.^\\{x)g\\2 + \\{x)Ay'fh M 
proving ()F.6p . The proof of ()F.7p is very similar. □ 



Step 2: Splitting of the integral. To bound KnAtp for n > 1, choose a smooth cut-off function x 
equal to 1 on i?(0, 1) and outside of i?(0, 2), then split the integral 

(F.8a) KnAi; =A j A,p{y)\^^^^^^^^x{x - y) dy 

(F-8b) +Aj ^^(y)L?^M_|gp(l - x{x - y)) dy . 

Observe that since ||xyl/||2 < ||xV/||2 + II/II2 we can replace yl / by V / in our estimate. 



GLOBAL SOLUTIONS FOR THE GRAVITY WATER WAVES EQUATION IN DIMENSION 3 41 

Step 3: Estimate of liF.8h\) : The case |x — y| > 1. In this region, the cancellation contained in 
h{x) — h{y) is not needed thus it suffices to replace {h{x) — h{y))'^^ by the general term given by 
the expansion of this power, namely h{xYh(y)'^"'~^ . We get 



For i = the above integral is given by 



(1 - xix - y))- 



dy. 



\h{y)\^^Ai,{y)V^ ^ ^%J^ dy = r* {A^ph^^, 

where F satisfies ||(x)-r||/^i < 1. Therefore 

x{r * (/IV'/i^")) = {xF) * {A^h^"") +^* {xA^h^'') 

Li L2 Li ^2 



L2 L2 

and the conclusion follows since 

\\xAijh^^h<\\xhh\\An^\\hf^-\ 

For £ > 1 the most problematic term is when V^^. hits h{x), since otherwise we have added decay 
in X — y. It reads 

A^{y){l - x{x - ^ Vh{x)h{xY-'n * {A^h^-') 



|L2 



\x — i/|2"+l 

where Fi decays at least like For ^ = 1 we have 



\x{Vh)iFi * (^V/i'"~')llL2 < ||V/i|U^|||x(ri * {A^ph^^-'] 



which can be estimated by 

x(ri * (yl^/i^"-!)) = (xFi) * + Fi * {xAiPh^"-^) 

L4/3 L4/3 LI L2 



L2 L2 

and since Lp'{{x)'^dx) ^ L^{dx), this yields the desired estimate as in the case £ = 0. For £ > 2 
the bound follows in an even simpler way, thus we skip this case. 

Step 4- Estimate for W.8a\) : The case \x — y\ ^ 1. We use Taylor's formula 

Mfl^ = VMx)f^ + V'h{x + t{y - - t)dt 

\x-y\ \y-x\ Jo \y-x\ 

to deal with the singular integral. More precisely, we substitute the above expression for each factor 
h{x) — h{y) in ()F.8aP : subsequently expanding the product gives 2^" terms. We start with the one 

involving ^V/i • ||zf| ) > writing it 

(F.9a) V,. I ^V(y)^^^^^ {vh{x) ■ dy 



\x — y\ \ \y — x\ 

(F.9b) =V, f^A^i^y)-A^{x))^^^(vh{x)-^^—^y^dy 

J \x-y\ V \y-x\j 

(F.9C) + V. {Ai^ix) I (Vhix) ■ ^) dy 
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The term ()F.9cP can be written as 

V ^p^2 \x-y\ \y-x\ \y-x\ ) 

= V{A^{x) ^ -fN^...N2^dNih- ■ ■ dN„h) 
which is easy to bound by (1F.6|) . To bound part (IF.9bj) . we note that if hits 71-0(2;) one gets 



VAi,{x) I 4^-^(vh{x) ■ ^^^y dy 

J \x-y\ V \y-x\J 

which is similar to part ()F.9cp . If V^; hits V/i(x) one gets 

{Ai,{y) - Ai,{x))4^^ ( V^h{x)-^^^^] ( Vh{x) ■ ' dy 

\x-y\ \ \y -x\J \ \y -x\J 

which can be written as a sum of terms of the type 

CV^h{x)Vh{xf"-'^A^ or V'^ h{x)V hixf""'^ K * Aip K an integrable kernel 

and thus can be treated by Lemma lF.31 

If V X hits the function of {x — y) we turn V^,. into — Vy and integrate by parts to get 



\y — x\ \ \y — x\ 

which can be written as a sum of terms of the type 

\/h[x)'^^K * Aifj K an integrable kernel 

and can be bounded by lemma IF.3[ 

All the remaining terms in the expansion of ()F.8aP after substituting the Taylor expansion involve 
at least one power of Jq V^/i . . . , i.e. 

V, [ Ai,{y)x{x-y)\x-y\'-^( [\^h{x + t{y - x)){l - t)dt] Vh{xf^-'dy 



for £ > 1 (here we abused notations by writing |x— instead of the correct multilinear expression 
in X — y involving scalar products). Regardless of what term is hit by V^^ we estimate the above in 
the following manner: 1) Split tp into low and high frequencies 

def 

^ = ''Plow + ^high = P>o^ + P>oi'; 

2) if ip is low frequency, estimate directly after putting the weight on -0; 3) if is high frequency, 
remove one derivative from -0 by integration by parts (thus getting Rip, R Riesz transform), then 
estimate after putting the weight on ip. Thus if ip is low frequency, we take the derivatives of h in 
L°° and get the bound 

\\{x)Ai,,M\hf^.,^ 

and if tp is high frequency, integrate by parts in y using Aip = VyRip and get the bound 

||(x)i?0high||2||/i||^3,oo 

The estimate follows since 

|Kx)Mow||2 + ||(x)i?0high||2 < \\{x)A'/^i;\\2 
This completes the proof of ()F.5ap . 
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Proof of (jF.5b|) in Claim IF. 21 The proof of the inequahty ()F.5bp is very similar to that of 
inequality (lF.5a[) . which was the object of the previous section. One additional ingredient is however 
needed, which reduces to the boundedness of the operator 

T-.Z^vj Z{y) - y) dy on L\ 

In order to prove it, observe that T is a singular integral operator if h belongs to W'^'°°. By the 
David and Journe Tl theorem, its boundedness over will follow from the belonging of T(l) and 
r*(l) to BMO. Since T is antisymmetric, it suffices to check that Tl belongs to BMO. But Tl 
is given by 

^ J l -y\2nll Xix-y)dy, 

which, as one sees easily, is bounded in provided h £ 1^^'°°. 

Remark. Equations (j3.ip can also be easily derive from ()F.3p and (|F.4p . Using the fact that 
A = '^Ridi, where Ri 's are the Riez potentials, one can show that 

thus from (jEll) p = A^ - ^[h'^A'^ip - 2hA{hAip) + A{h'^ATp)] + and from (IF^ 

G{h)i; = yl^ - V • (HViIj) - A{hAi;) - ]^[A{h^A^^) + A^{h^Ai)) - 2A{hA{hAil:))] + . . . , 



which is the first part of (13. ip . The second part follows immediately from the first 
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